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EXECUTIVE  SUMMARY 


The  Tomahawk  Land- Attack  Missile  (TLAM)  is  the  long-range  accuracy  strategic 
weapon  of  choice  in  strike  warfare  against  strategic  targets  for  United  States  (U.S.) 
military  forces.  When  the  U.S.  National  Command  Authority,  formed  by  the  President 
and  the  Secretary  of  Defense,  authorizes  a  TLAM  strike,  the  authorization  passes  down 
through  the  chain  of  command  via  the  Chief  of  Naval  Operations,  a  regional  fleet 
Commander-in-Chief,  and  ultimately  to  the  Battle  Group  Tomahawk  Strike  Coordinator 
(TSC). 

Predesignation  is  the  process  of  determining  which  ship  or  submarine  will  fire 
which  TLAM  missiles  in  support  of  an  authorized  attack  upon  specified  targets. 
Predesignation  has  two  phases.  In  Phase  1,  the  TSC  allocates  targets  to  ships  and 
submarines.  In  Phase  2,  tasks  are  allocated  to  specific  TLAMs  onboard  individual  ships 
and  submarines.  Currently,  the  TLAM  is  mostly  being  used  to  attack  a  fixed  stationary 
target  such  as  a  building,  but  under  the  conditions  of  modem  warfare,  the  TLAM  must  be 
able  to  attack  a  mobile  target  on  short  notice.  Our  goal  is  to  provide  the  TSC  with  an 
automatic  predesignation  procedure  that  enables  tasks  to  be  assigned  quickly  and 
efficiently. 

Over  the  last  five  years,  faculty  and  students  at  the  Naval  Postgraduate  School 
(NPS)  Operations  Research  Department  have  addressed  this  problem  under  the 
sponsorship  of  the  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD).  A 
series  of  student  theses  and  other  reports  document  this  effort.  As  a  result,  a  multi¬ 
objective  (MO)  mixed- integer  programming  (MIP)  model  was  developed  by  the  end  of 
1999.  This  model  represented  an  important  challenge  in  terms  of  solvability.  In  fact, 
previous  work  determined  that  real-size  problems  cannot  be  solved  within  reasonable 
optimality  tolerance  in  affordable  time.  This  difficulty  was  addressed  in  subsequent 
years  by  means  of  fast  heuristic  algorithms  at  the  expense  of  solution  quality  and,  to  some 
extent,  flexibility  in  the  set  up  of  goal  priorities  and  other  potential  redesigns  in  the 
problem  specifications.  As  the  NPS  team’s  knowledge  of  TLAM  predesignation 

xxiii 


evolved,  the  heuristies  addressed  features  of  the  problem  that  were  not  eonsidered  in  the 
earliest  optimization  models.  The  latest  heuristie  algorithm  has  been  reeently 
implemented  by  professional  programmers  at  NSWCDD,  adhering  to  algorithm 
speeifioations  by  NPS  faeulty.  This  thesis  bridges  the  gap  between  the  latest 
speeifieations  implemented  in  the  heuristie  algorithms  and  the  outdated  mathematieal 
formulation.  We  show  that  the  revised  formulation  ean  be  solved  in  reasonable  time  for 
the  existing  test  eases,  ineluding  those  that  ineorporate  reeent  features  that  were  not 
ineorporated  in  earlier  versions  of  the  formulation. 

Speeifie  aeeomplishments  of  this  thesis  are: 

•  Consolidation  of  all  problem  speeifieations  to  date. 

•  Ineorporation  of  all  speeifioations  into  a  revised  MO  MIP  formulation.  In  partioular, 
the  following  oonsiderations  have  been  added  to  the  optimization: 

•  Submarines  and  their  oharaoteristios:  torpedo  tubes,  missiles  in  the 
torpedo  room,  oanister  oapsules,  eto. 

•  Target  oonoept  as  a  set  of  multiple  tasks  sharing  oertain  joint 
oonstraints. 

•  Multiple  missiles  might  be  required  for  every  task-part. 

•  TSC’s  manual  allooation  missiles  to  speeifie  oells,  torpedo  tubes  or 
oanister  oapsules. 

•  TSC  speoifioation  regarding  whether  multiple  missiles  for  a  given 
task-part  must  be  assigned  to  the  same  firing  unit,  or  separating 
missiles  aoross  multiple  firing  units  is  allowed. 

•  TSC  speoifioation  regarding  whether  baok-up  task-parts  ean  be 
allooated  to  submarines  or  not. 

•  Ineorporation  of  ghost  task-parts  and  their  oharaoteristios. 

•  Exploration  of  alternative  teohniques  to  obtain  aoourate  bounds  on  large-size  test 
oases.  We  oompare  the  results  against  those  provided  by  the  heuristie  to  provide 
empirio  validation  for  using  the  heuristie  as  a  reliable  tool  to  support  fast  TEAM 
predesignation. 
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I.  INTRODUCTION 


A.  BACKGROUND 

1.  History 

After  the  signing  of  first  Strategic  Arms  Limitation  Treaty,  the  strategic  cruise 
missile  was  born.  Until  January  6,  1977,  the  United  States  (U.S.)  Navy  and  U.S.  Air 
Force  maintained  their  own  projects  independently.  On  that  date,  the  Defense  Systems 
Acquisition  Review  Council  authorized  a  full-scale  engineering  and  testing  program  for 
air-launched,  anti-surface,  theater,  and  nuclear  versions  of  the  cruise  missile.  The  Joint 
Cruise  Missiles  Project  Office  was  formed  to  manage  cruise  missile  development,  and  it 
became  a  joint  U.S.  Navy  and  Air  Force  project.  Today,  the  U.S.  Navy  continues  as  a 
lead  agency  in  cruise  missile  program  management.  [Reynolds,  1998.] 

The  first  Tomahawk  Land  Attack  Missile  (TLAM)  was  deployed  on  board  a  Los 
Angeles  class  fast-attack  submarine  in  November  1983.  Its  first  operational  use  was 
during  Operation  Desert  Storm  in  1991.  Since  its  first  combat  use,  it  has  become  the 
weapon  of  choice  in  strike  warfare  for  U.S.  Military  forces.  The  governments  of  the  U.S. 
and  United  Kingdom  (U.K.)  signed  a  Foreign  Military  Sales  Agreement  for  the 
acquisition  of  65  missiles  in  1995,  recorded  as  the  first  sale  of  TLAMs  to  a  foreign 
country  [Navy  Office  of  Information,  2004.].  After  a  successful  launch  and  live  warhead 
test  in  November  1998,  the  U.K.  declared  operational  capability.  TLAMs  can  be 
currently  launched  from  U.S.  Navy  surface  ships  and  from  U.S.  Navy  and  Royal  Navy 
submarines. 

2,  Characteristics  of  Tomahawk  Land  Attack  Missile 

The  TLAM  (Figure  1)  is  a  long  range,  subsonic  cruise  missile  designed  for  land 
attack  warfare.  It  was  designed  to  hit  strategic  targets  day  and  night  with  a  minimum  of 
collateral  damage  by  flying  extremely  low  altitudes  at  high  subsonic  speeds  and  piloting 
over  an  evasive  route  by  several  tailored  mission-guidance  systems.  It  is  able  to  convey  a 
1,000  pound  warhead  to  a  range  of  approximately  l,000nm  with  pinpoint  accuracy. 
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Figure  1.  Tomahawk  Land  Attack  Missile  [Onwardoverland,  2004] 


Some  TLAM  characteristics  [Navy  Office  of  Information,  2004.]  are  as  follows: 
Unit  Cost:  Approximately  $569,000  (FY99  $) 

Length:  18  feet  3  inches  (5.56  m);  with  booster:  20  feet  6  inches  (6.25  m) 

Weight:  2,900  pounds  (1,315.44  kg);  3,500  pounds  (1,587.6  kg)  with  booster 
Diameter:  20.4  inches  (51.81  cm) 

Wing  Span:  8  feet  9  inches  (2.67  m) 

Speed:  Subsonic  -  about  550  mph  (880  km/h) 

Power  Plant:  Block  II/III  TLAM- A,  C,  D:  Williams  International  FI 07  cruise  Block 
IV  TLAM-E:  Williams  International  F415  cruise  turbo-jet  engine,  ARC  solid- fuel 
booster 

Range:  Block  II  TLAM-A:  1,350  nm  (1,500  statute  miles,  2,500  km) 

•  Block  III  TLAM-C:900  nm  (1,000  statute  miles,  1,600  km) 

•  Block  III  TLAM-D:700  nm  (800  statute  miles,  1,250  km) 

•  Block  IV  TLAM-E:900  nm  (1,000  statute  miles,  1,600  km) 
Guidance  System: 

•  Block  II  TEAM-A:  Inertial  Navigation  System  (INS)  Terrain 
Contour  Matching  (TERCOM) 

•  Block  III  TLAM-C,  D:INS  TERCOM  Digital  Scene  Matching 
Area  Correlation  Global  Positioning  Satellite  System  (GPS) 

•  Block  III  TEAM-C,  D:INS  TERCOM  GPS 

Warheads: 

•  Block  II  TEAM-N:  W80  nuclear  warhead 


•  Block  III  TLAM-C:  1,000-pound  class  unitary  warhead 

•  Block  IV  TLAM-E:  1 ,000-pound  class  unitary  warhead 

•  Block  III  TLAM-D:  Conventional  submunitions  dispenser  with 
combined  effect  bomblets 


•  Date  Deployed: 

•  Bloek  II  TEAM -A:  Initial  Operational  Capability  (IOC)  in  1984 

•  Bloek  III:  IOC  in  1994 


•  Bloek  IV:  Expected  2004 

B,  FIRING  UNITS 


TEAMS  ean  be  fired  from  ships  and  submarines.  Some  charaeteristies  of  these 
firing  units  are  shown  in  Eigure  2. 
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Eigure  2.  TEAM  Eiring  Units. 

Ship  Vertieal  Eaunehing  System  (VES)  offer  more  storage  capaeity  than  submarines, 
which  may  or  may  not  have  a  Canister  Eaunehing  System  (CES).  [Arnold,  2000] 

1,  Ships 

Ships  can  store  more  TEAMs  than  submarines.  They  can  also  maintain  better 
operational  communication  capabilities  in  order  to  eoordinate  TEAM  missions  with  other 
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firing  units  and  the  Tomahawk  Strike  Coordinator  (TSC).  Ships  can  be  considered  as  the 
main  TLAM  firing  units.  [Arnold,  2000.] 

TLAMs  on  ships  are  launched  from  a  Vertical  Launching  System  (VLS). 
Depending  on  its  type,  a  ship  may  have  a  forward  launcher  only  or  both  aft  and  forward 
launchers.  A  launcher  can  be  half-size  or  full-size  depending  on  the  type  of  ship.  The 
basic  component  of  a  launcher  is  called  a  module.  A  module  provides  power  and 
computer  connections  for  each  weapon  to  be  fired  from  the  cells  in  the  module.  Modules 
may  be  of  standard  type  or  crane  type.  The  standard  module  is  divided  into  two  half¬ 
modules  and  each  half-module  is  divided  into  four  cells.  In  the  crane  module,  the  crane 
occupies  cells  “6”  to  “8”  (see  Figure  3).  Thus,  the  crane  module  has  a  total  of  two  half¬ 
modules  and  five  cells.  Each  cell  can  store  only  one  TLAM.  Half-size  launchers  have 
four  modules,  eight  half-modules,  29  cells  and  can  store  up  to  29  TLAMs,  whereas  full- 
size  launchers  have  eight  modules,  16  half-modules,  61  cells  and  can  store  up  to  61 
TLAMs  (see  Figure  4).  [Arnold,  2000.] 


Figure  3.  Standard  and  Crane  Modules. 

The  basic  unit  of  the  VFS  is  a  module.  The  standard  module  is  divided  into  two  half¬ 
modules  consisting  of  cells  “1-4”  and  “5-8”  respectively.  In  a  crane  module,  cells  “6-8” 
are  occupied  by  the  crane  and  cell  “5”  is  the  only  one  available  in  the  half-module,  while 
the  other  half-module  consists  of  cells  “1-4.”  [Arnold,  2000] 
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Aft  Launcher 


Forward  Launcher 


VERTICAL  LAUNCHING  SYSTEM  (FULL)  VERTICAL  LAUNCHING  SYSTEM  (HALF) 

Figure  4.  Full-Size  and  Flalf-Size  VLS  Layouts. 

A  full  launeher  eonsists  of  eight  modules  whereas  a  half-size  launeher  eonsists  of  four 
modules.  Both  launehers  have  one  crane  module  with  5  cells.  The  circled  cells  are 
designated  as  “Aft,  3,  5,”  and  “Forward,  8,  4.”  [Arnold,  2000] 


Salvo  size  refers  to  the  number  of  TLAMs  that  can  be  aligned  and  fired  by  a  firing 
unit  during  any  conflict  time  window  (typically  a  45-minute  period).  Due  to  VLS’s 
design  and  power  limitations,  only  one  missile  can  be  fired  from  a  half-module  during  a 
window.  The  tasks  that  are  to  be  aligned  within  the  same  window  are  said  to  be  near- 
simultaneous  or  to  conflict  with  each  other  (see  Section  D  for  details).  Because  of  half¬ 
module  power  constraints,  a  full-size  launcher  has  a  maximum  salvo  capability  of  16 
missiles  during  any  window,  while  a  half-size  launcher’s  salvo  is  eight  missiles.  [Arnold, 
2000.] 
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2,  Submarines 

Submarines  can  approach  targets  stealthily,  but  they  must  be  close  to  the  surface 
with  antennas  above  the  water  when  firing  TLAMs.  They  reveal  their  location  after  the 
first  firing.  Submarines  may  have  communication  problems  with  the  TSC  in  the 
operational  area:  for  example,  a  TSC’s  short-notice  order  may  not  be  received  on  time 
due  to  limited  communication  capability  for  security  reasons.  [Arnold,  2000.] 

Submarines  have  two  assets  to  align  TLAMs:  torpedo  tubes  and  canister  capsules. 
All  U.S.  Navy  attack  submarines  (SSNs)  have  the  capability  of  firing  TLAMs  from  their 
torpedo  tubes  [See  Figure  5].  A  TLAM  can  be  loaded  into  a  torpedo  tube,  and  then 
aligned  and  fired  like  a  torpedo.  Torpedo  tubes  can  be  reloaded  with  TLAMs  carried  on 
the  submarine.  All  torpedo  tubes  on  a  submarine  are  capable  of  launching  TLAMs. 
However,  in  order  to  maintain  enough  self-defense  capability,  a  commanding  officer 
might  prefer  not  to  dedicate  all  torpedo  tubes  to  TLAM  missions.  [Arnold,  2000.] 

In  addition  to  torpedo  tubes,  Los  Angeles  class  submarines  with  hull  number  of 
719  (or  higher)  have  a  Capsule  Launching  System  (CLS)  that  provides  12  canister 
capsules  exclusively  for  firing  TLAMs  (see  Figure  6)  [Wingeart,  2001].  The  CLS  on  a 
submarine  is  similar  to  the  surface  ship’s  VLS  and  cannot  be  reloaded  at  sea.  Unlike  a 
ship  VLS,  a  submarine  CLS  is  not  restricted  by  power  constraints,  so  all  missiles  in  CLS 
cells  may  be  aligned,  if  necessary,  within  a  conflict  time  window. 
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Figure  5.  Cutaway  View  of  a  Torpedo  Tube  Room. 

This  Los  Angeles  Class  Submarine  has  four  torpedo  tubes  (TB1,...,TB4).  The  ordnance 
is  stored  in  the  torpedo  room.  Five  CII  TLAMs  and  ten  CIII  TLAMs  ordnance  are  stored 
in  this  torpedo  room.  [Arnold,  2000] 


Figure  6.  Capsule  Launching  System  . 

This  CLS  has  12  capsules;  each  can  launch  all  current  blocks  and  variants  of  TLAMs  and 
is  not  limited  by  any  power  constraint.  [Arnold,  2000] 

C.  PREDESIGNATION 

Predesignation  is  the  assignment  of  land-attack  missiles  (on  ships  and 
submarines)  to  target  aim-points.  This  decision  process  has  two  phases:  (1)  the  allocation 
of  the  TLAMs  to  the  firing  units,  considering  all  firing  units  and  tasks  simultaneously. 
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and  (2)  given  the  firing  unit  alloeation  from  Phase  1,  a  speeifie  alloeation  of  tasks  to 
missiles.  The  models  presented  in  earlier  works  and  in  this  thesis  address  the  final 
alloeation  to  specific  missiles  on  firing  units,  which  is  the  ultimate  goal  of  Phase  2. 

When  the  U.S.  National  Command  Authority  (NCA),  formed  by  U.S.  President 
and  the  Secretary  of  Defense,  authorizes  a  strike,  the  authorization  passes  down  through 
the  chain  of  command  via  the  Chief  of  Naval  Operations  (CNO)  and  a  regional  fleet 
Commander-in-Chief  (CINC)  and  eventually  to  the  Battle  Group  TSC.  To  begin  Phase  1, 
the  NCA  or  CINC  provides  the  TSC  the  following  information;  a  prioritized  list  of 
targets,  a  list  of  aim-points  at  the  target,  the  desirable  and  feasible  TLAM  types,  the 
number  of  missiles  per  aim-points,  missile  time  on  the  target,  and  missile  assignment 
redundancy  requirement.  [Arnold,  2000.] 

The  TSC  predesignates  the  aim-points  to  firing  units  by  selecting  a  set  of  pre¬ 
planned  missions  that  meet  NCA  or  CINC  criteria.  A  pre-planned  mission  defines  the 
target,  aim-points,  block,  warhead  and  launch  area  [Arnold,  2000].  Predesignation 
considers  the  geo-graphical  nearness  of  the  candidate  firing  units  to  aim-points,  the 
inventory  and  location  of  TLAMs  aboard  each  firing  unit,  the  limitations  of  preparing  and 
firing  particular  missiles,  the  flight  route  coordination  among  TLAMs  and  other  tactical 
concerns  [Brown  et  ah,  2001].  The  other  important  issue  to  be  considered  by  the  TSC  is 
to  leave  the  combat  units  with  a  maximum  residual  salvo  size  after  a  strike,  individually 
and  as  a  battle  group.  The  TSC  tries  to  preserve  as  many  remaining  salvos  as  possible  for 
each  missile  type  on  each  firing  unit  that  will  remain  with  the  battle  group  in  the 
operation  area.  The  TSC  also  avoids  predesignations  that  interfere  with  other  duties  of 
the  firing  units  and  attempts  to  expend  TLAMs  from  the  firing  units  that  will  soon  be 
departing  the  theater. 

D,  STRIKE  PLANNING  CONSIDERATIONS 

1.  General  Considerations 

A  10-character  alphanumeric  code,  known  as  the  Basic  Encyclopedia  Number 
(BEN),  is  used  to  designate  each  target  chosen  for  attack.  [Wingeart,  2001.]  Eor 
example,  the  BEN  of  “Target  1”  may  be  3004005000.  The  TSC  determines  the  aim- 
points  to  refine  the  impact  points  of  TLAMs  on  the  target  and  the  best  TLAMs  suited  to 
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attack  the  targets.  In  addition,  the  TSC  must  determine  the  missions  required  to  hit  the 
aim-points  he  or  she  designates,  the  variants  and  total  number  of  TLAMs  to  fire,  and  the 
desired  exeeution  time  at  each  aim-point. 

An  1 1 -digit  number  known  as  the  Mission  Identification  Number  (MIN)  speeifies 
eaeh  task  (also  called  mission)  and  associated  list  of  missile  variants  that  can  be  used  to 
aceomplish  the  assigned  tasks.  [Arnold,  2000.]  Feasible  missile  variants  for  the  task  are 
listed  on  the  Missile  Mission  Matching  (M  )  list.  Every  task  has  a  list  that  indieates  the 
relative  desirability  of  all  capable  missile  variants.  The  TSC  should  use  the  least  capable 
missile  in  order  to  maximize  residual  salvo  capabilities,  as  long  as  this  ehoice  does  not 
interfere  with  other  priorities.  For  instance,  if  the  deeision  should  be  either  to  select  a 
Block  C-II  or  a  Block  C-III  TEAM,  a  Block  C-II  is  preferred  sinee  it  has  a  less-capable 
warhead  and  shorter  range  than  the  Bloek  C-III  missile. 

2,  Near-Simultaneous  Tasks 

As  explained  in  Section  B,  the  number  of  TLAMs  fired  from  a  VLS  launcher  on  a 
ship  is  restricted  by  a  power  constraint,  which  allows  only  one  TEAM  to  be  aligned  at 
once  per  half-module.  From  torpedo  tubes  on  submarines  this  number  is  precisely  the 
number  of  torpedo  tubes.  In  this  context,  at  once  means  within  a  pre-specified 
eonflicting  time  window,  typieally  45  minutes. 

A  task  has  a  launch  time  window  defined  by  an  earliest  time  to  launeh  (ETL)  and 
latest  time  to  launch  (ETL).  [Kubu,  2001.]  These  are  derived  from  the  required  arrival 
time  at  the  aim-point  and  missile  flight  parameters.  Although  it  would  be  desirable  to 
fully  employ  this  information  to  add  flexibility  to  the  predesignation  process,  our 
approaeh  is  a  simplification  of  the  time  window  specification.  This  is  due  to  the  fact  that 
we  will  not  optimize  the  exact  time  within  the  ETL-LTL  window.  Instead,  we  eonsider 
an  average  time  to  launch  (ATE),  calculated  as  the  average  of  earliest  and  latest  times  to 
launeh.  If  the  average  time  to  launch  between  any  two  tasks  is  less  than  the  given 
eonflieting  time  window,  they  are  eonsidered  near-simultaneous  tasks,  and  they  are 
incorporated  into  a  so-ealled  conflict  task  set  or  near-simultaneous  task  set.  A  task  may 
be  in  more  than  one  of  these  sets. 
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The  following  example  elarifies  the  eoneept  of  oonfliet  sets  and  near- 
simultaneous  tasks  (see  Figure  7).  In  this  example,  we  assume  the  eonflieting  time 
window  is  45  minutes.  While  some  tasks  are  in  only  one  near-simultaneous  task  set, 
others  are  in  more  than  one  of  these  sets,  depending  on  their  ATLs.  Assuming  that  C-II 
and  C-III  are  valid  weapons  for  all  of  our  tasks,  by  eombining  Figures  7  and  8  we 
eonelude  that,  for  example,  the  TSC  eannot  alloeate  tasks  “150-150-15000”  and  “155- 
155-15500”  to  the  module  in  Figure  8  due  to  power  eonstraints.  However,  the  TSC  may 
alloeate  tasks  “155-155-15500,”  “160-160-16000”  and  “175-175-17500”  to  this  module. 


TASK  CODE 
(MIN) 

ETL  (Earliest 
Time  to 
Launch) 

LTL  (Latest 
Time  to 
Launch) 

ATL  (Average 
Time  to 
Launch) 

Average 
Time  (min) 
Since 
8/19/2000 
00:00:00  AM 

Conflict 
Task  Set 
Number  (s) 

150-150-15000 

08/19/2000  00:00 

08/19/2000  00:10 

08/19/2000  00:05 

5 

1 

155-155-15500 

08/19/2000  00:00 

08/19/2000  00:10 

08/19/2000  00:05 

5 

1 

160-160-16000 

08/19/2000  00:30 

08/19/2000  00:40 

08/19/2000  00:35 

35 

1,2 

165-165-16500 

08/19/2000  00:30 

08/19/2000  00:40 

08/19/2000  00:35 

35 

1,2 

170-170-17000 

08/19/2000  00:30 

08/19/2000  00:40 

08/19/2000  00:35 

35 

1,2 

175-175-17500 

08/19/2000  01:00 

08/19/2000  01:10 

08/19/2000  01:05 

65 

2,3 

180-180-18000 

08/19/2000  01:00 

08/19/2000  01:10 

08/19/2000  01:05 

65 

2,3 

185-185-18500 

08/19/2000  01:00 

08/19/2000  01:10 

08/19/2000  01:05 

65 

2,3 

190-190-17000 

08/19/2000  01:00 

08/19/2000  01:10 

08/19/2000  01:05 

65 

2,3 

225-225-22500 

08/19/2000  01:30 

08/19/2000  01:40 

08/19/2000  01:35 

95 

3 

230-230-23000 

08/19/2000  01:30 

08/19/2000  01:40 

08/19/2000  01:35 

95 

3 

235-235-23500 

08/19/2000  01:30 

08/19/2000  01:40 

08/19/2000  01:35 

95 

3 

Figure  7.  Confheting  Sets  of  Tasks. 

ATL  is  ealeulated  by  averaging  ETL  and  LTL.  We  assume  that  two  tasks  whose  ATLs 
are  within  45  minutes  of  eaeh  other  are  near-simultaneous.  In  this  example,  tasks  “150- 
150-15000”  through  “170-170-17000”  are  all  eonflieting  with  eaeh  other.  Thus,  they 
belong  to  oonfliet  set  #1.  The  same  ooours  for  tasks  “160-160-16000”  through  “190-190- 
17000,”  which  deflnes  conflict  set  #2.  Notice  that  some  tasks  belong  to  more  than  one 
conflict  set  [Excerpt  from  TEAM  Documentation,  2002]. 
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Figure  8.  Module  with  an  Empty  Half-Module  and  Full  Half-Module. 

The  two  C-IIs  and  two  C-IIIs  in  the  upper  half-module  ean  be  used  for  non-near- 
simultaneous  tasks.  For  example,  only  one  of  tasks  “150-150-15000”  and  “155-155- 
15500”  (from  Figure  7)  ean  be  alloeated  to  the  upper  half-module. 

3.  Missile  Redundancy 

One  of  the  most  important  eonsiderations  for  the  TSC  is  to  deeide  the  amount  of 
redundaney  for  eaeh  task  of  a  target.  Intuitively,  the  TSC  wants  to  alloeate  more  missiles 
to  high  priority  targets  in  order  to  ensure  their  destruetion  in  ease  some  missiles  fail  to 
launeh  or  miss  the  target.  This  leads  to  eonsideration  of  primary,  ready-spare,  and  baek- 
up  task-parts,  as  well  as  a  speeial  type  of  task-parts  ealled  ghost.  A  primary  task-part  is 
designated  to  be  launehed  at  a  target,  whereas  ready-spare  and  baek-up  task-parts  are 
only  exeeuted  in  the  ease  of  primary  task-part  failure.  Primary  and  ready-spare  task-parts 
must  be  assigned  to  the  same  firing  unit,  whereas  baek-up  task-parts  must  be  assigned  to 
a  different  firing  unit  than  the  primary  and  ready-spare  (if  any)  task-parts.  [Arnold,  2000.] 

Multiple  missiles  may  be  required  for  eaeh  task-part.  The  ready-spare  and  baek- 
up  task-parts  are  aligned  simultaneously  with  the  primary  task-part.  Ready-spare  task- 
parts  are  launehed  immediately  if  the  missile  assoeiated  with  the  primary  task-part  fails, 
and  baek-up  task-parts  are  launehed  if  the  missiles  assoeiated  with  the  primary  task-part 
and  ready-spare  (if  any)  fail.  Ghost  task-parts  ean  only  be  assigned  to  a  speeial  type  of 
tasks,  ealled  ghost  tasks.  A  ghost  task  enables  one  missile  to  provide  redundaney  for 
multiple  tasks  and  multiple  task-parts  eoneurrently.  A  ghost  task  is  fired  when  ordered 
by  the  TSC.  For  some  ealeulations  (sueh  as  expeeted  residual  salvo),  we  assume  that  all 
missiles  assoeiated  with  primary  task-parts  are  fired,  and  those  assoeiated  with  ready- 
spare,  baek-up  and  ghost  task-parts  are  not.  [Arnold,  2000.] 
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4,  Other  Considerations 

There  are  some  doctrinal  constraints  applying  to  targets  and  tasks.  Each  target  is 
designated  either  as  restricted  or  unrestricted.  All  the  task-parts  in  tasks  associated  with 
unrestricted  targets  may  be  spread  among  the  geo-feasible  firing  units.  Restricted  targets 
have  their  tasks  and  associated  task-parts  allocated  to  a  limited  number  of  firing  units, 
which  is  specified  by  the  TSC.  In  addition,  if  any  task-part  of  a  task  in  a  restricted  target 
is  left  unassigned,  all  the  other  task-parts  of  that  task  will  also  be  considered  as 
unassigned.  [Arnold,  2000.] 

For  each  task-part  with  multiple  missile  requirements  the  TSC  specifies  whether 
the  missiles  for  the  task-part  are  separable  (formerly  referred  to  as  split)  or  not  across 
firing  units.  If  missiles  of  any  task-part  are  non-separable,  all  missiles  of  that  task-part 
should  be  assigned  to  the  same  firing  unit.  For  example,  in  Figure  9,  “Target  A”  has  two 
tasks:  “Task  1”  and  “Task  2.”  These  tasks  have  their  own  task-parts  (e.g.,  there  are  two 
primary  missiles  required  for  “Task  1”). 


TARGET  A 


TASK  1 


TASK  2 


f 


1 


Primary 

Ready- spare 

Back-up 

Primary 

Ready- spare 

Back-up 

(SEPARABLE) 

1  1 

i 

i 

(NON-SEPARABLE) 

1  1 

i  1 

1  i 

PRlI  |PR2| 

RSI 

BUI 

1  PRl  PR2| 

RSI  RS2 

BUI  BU2 

Allocated  to  Allocated  to 
Firing  Unit  1  Firing  Unit  2 


Allocated  to 
Firing  Unit  1 


Figure  9.  Separable  and  Non-separable  Tasking. 

In  the  figure,  primary  task-parts  in  “Task  I”  are  deemed  “separable,”  but  primary  task- 
parts  in  “Task  2”  are  “non-separable.”  Thus,  “PRl”  and  “PR2”  in  “Task  1”  can  be 
allocated  either  to  the  same  firing  unit  or  to  different  firing  units,  whereas  “PRl”  and 
“PR2”  in  “Task  2”  must  be  allocated  to  the  same  firing  unit. 
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If  the  primary  task-part  in  “Task  1”  is  deemed  separable,  we  ean  alloeate  both 
“PRl”  and  “PR2”  to  either  the  same  firing  unit  or  to  different  firing  units.  On  the  other 
hand,  if  the  primary  task-part  in  “Task  2”  is  deemed  non-separable,  both  “PRl”  and 
“PR2”  must  be  alloeated  to  the  same  firing  unit. 

E.  SUMMARY  OF  PAST  WORK 

Sinee  1998,  a  group  of  faculty  and  students  at  the  Naval  Postgraduate  School 
(NPS)  Operations  Research  Department  has  contributed  to  the  development  of  a  decision 
support  algorithm  for  Tomahawk  Predesignation,  under  the  sponsorship  of  the  Naval 
Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD).  A  series  of  student  theses 
documents  this  effort.  Most  of  these  theses  were  advised  by  Professor  Alexandra  M. 
Newman,  now  with  the  Colorado  School  of  Mines. 

LT  Scott  Kuykendall  [1998]  breaks  ground  by  addressing  Tomahawk 
predesignation  with  a  mixed-integer  programming  (MIP)  formulation  with  elegant 
simplicity.  LT  Kuykendall’s  fleet  predesignation  experience  motivated  him  as  follows: 
(i)  current  paper-and-pencil  predesignation  process  is  inefficient  and  error-prone,  and  (ii) 
current  predesignation  software  is  neither  user-friendly  nor  satisfactory,  since  it  produces 
predesignations  that  can  be  trivially  improved  by  inspection  [Wu  et  ah,  2002]. 
Kuykendall’s  work  considers  a  single-ship  scenario  and  aims  to  maximize  residual  strike 
capability  (also  referred  to  as  follow-on  salvo  or  residual  salvo).  He  proves  this  scheme 
yields  better  solutions  than  the  manual  predesignation  method;  yet,  his  MIP  does  not 
address  all  relevant  objectives. 

LT  Brian  Kirk  [1999]  develops  a  multi-objective  (MO)  MIP  for  automatic  surface 
ship  TLAM  Phase  I  predesignation  at  the  fleet  level.  His  model  is  based  on  eight 
objectives,  reviewed  by  Tactical  Training  Groups  Atlantic  and  Pacific,  presented  below 
in  descending  order  of  importance  [Wu  et  ah,  2002]: 

•  Minimize  the  number  of  unassigned  tasks  (actually  task-parts). 

•  Minimize  the  use  of  units  already  occupied  with  other  operations  or  not  in 
geo-graphical  proximity  to  the  strike. 

•  Maximize  the  allocation  of  tasks  across  so-called  expend  firing  units. 

•  Use  the  missiles  from  firing  units  designated  as  non-expend  so  as  to  level  the 
number  of  missiles  remaining  across  these  firing  units. 
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•  Spread  primary  task-part  missiles  aeross  as  many  tiring  units  as  possible  to 
prevent  single-point  failures  among  primary  task-part  alloeation. 

•  Spread  baek-up  task-part  missiles  aeross  as  many  firing  units  as  possible  to 
prevent  single-point  failures  among  back-up  task-part  allocation. 

•  Use  the  least  capable  missile  possible  for  each  mission. 

•  Maximize  the  residual  salvo  capability. 

Kirk  applies  three  solution  methods.  First,  he  attempts  to  solve  the  model  as  a 
single,  monolithic  problem  consisting  of  a  combination  all  the  objectives  into  a  single 
objective  function,  where  higher-priority  goals  have  greater  weights.  He  finds  that  this 
technique  does  not  produce  sensible  answers  for  large  scenarios,  and  solution  times  are 
prohibitive.  Next,  he  implements  a  Hierarchical  Restriction  by  solving  for  a  goal  at  a 
time.  He  achieves  solutions  of  reasonable  quality  by  using  this  technique  but  solution 
times  are  still  too  large  to  be  operationally  useful.  In  order  to  lessen  the  solution  times, 
Kirk  implements  a  Heuristic  Hierarchical  Restriction  method  on  the  MO  MIP,  which 
produces  faster  solutions  compared  to  those  of  previous  techniques,  but  not  fast  enough 
to  be  used  in  an  operational  setting.  [Wu  et  ah,  2002.] 

LT  Bertram  Hodge  [1999]  introduces  a  new  heuristic,  based  on  a  prioritized  target 
list  that  he  uses  to  mimic  the  optimal  decisions  of  Kirk’s  most  comprehensive  model. 
[Brown  et  ah,  2001.]  This  fast,  heuristic  algorithm  selects  firing  units  and  then 
predesignates  targets  from  the  list  in  a  single  pass.  Hodge  uses  an  algorithm  to  select 
ships  and  then  performs  task  allocations  based  on  a  single  pass  of  a  prioritized  task  list. 
His  heuristic  produces  solutions  roughly  20%  worse  than  Kirk’s  Hierarchical  Restriction 
for  all  objective  function  values.  As  a  result,  his  model  reduces  solution  times  at  the  cost 
of  solution  quality. 

LT  Andrew  Arnold  [2000]  improves  the  performance  of  Hodge’s  heuristic  by 
modifying  the  heuristic  algorithm  to  improve  solution  quality.  Arnold’s  heuristic  extends 
Hodge’s  fast  heuristic  to  include  submarines  as  firing  units  and  ghost  missions  as  task- 
parts.  Arnold  introduces  the  concepts  of  Target  (as  a  group  of  tasks  with  common 
constraints),  as  well  as  the  possibility  of  multiple  missiles  assigned  for  any  task-part.  He 
also  introduces  the  concept  of  missile  separability  by  task-parts,  and  allows  the  TSC  to 


specify  part  of  the  allocation.  [Wu  et  ah,  2002] 
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LT  Justin  A.  Kubu  [2001]  enhances  Arnold’s  heuristic  by  introducing  diagnoser 
and  prescriber  modules  that  provide  the  TSC  with  possible  reasons  why  a  task  cannot  be 
allocated  by  the  heuristic  algorithm.  Although  the  heuristic  cannot  determine  whether  a 
feasible  allocation  under  the  given  conditions  exists,  it  prescribes  modifications  to 
facilitate  a  complete  allocation  of  tasks.  Prescriptions  include,  for  example,  shifting  the 
launch  time  of  a  task  for  situations  in  which  the  heuristic  predesignation  does  not  provide 
a  complete  allocation. 

LT  Paul  H.  Wingeart  [2001]  enhances  Arnold’s  heuristic  by  implementing  two 
computer-assisted  allocation  methods:  (i)  task-part  to  firing  unit,  and  (ii)  firing  unit  to 
task-part.  Wingeart  validates  his  heuristic  solutions  with  exercise  data  TEXAS 
THUNDER  00-6  and  RIMPAC  00-3  and  shows  that  automated  allocation  yields 
solutions  superior  to  those  of  manual  predesignation  in  these  fleet  exercises. 

Consolidated  specifications  for  the  problem  data  and  heuristic  algorithm  can  be 
found  in  TEAM  documentation  [2002]. 

F.  THESIS  OUTLINE 

The  rest  of  the  thesis  is  organized  as  follows:  Chapter  11  introduces  the  problem 
specifications  that  serve  as  basis  for  our  mathematical  formulation.  Then  we  introduce 
the  problem  formulation  as  a  MIP  MO  optimization  model.  We  present  the  so-called 
hierarchical  and  combined-goal  formulations  for  the  problem.  Chapter  111  summarizes 
the  results  for  all  test  cases  from  our  implementation  in  GAMS  [Brooke  et  ah,  1996].  We 
also  describe  a  variable  partition  technique  called  Fix-and-Relax  to  accelerate 
convergence  in  some  cases.  Chapter  IV  presents  our  conclusions  and  recommendations 
for  future  work. 
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II.  MIXED-INTEGER  LINEAR  OPTIMIZATION  MODEL 


In  this  chapter,  we  present  the  mathematieal  formulation  of  our  TLAM 

predesignation  problem  as  a  MO  MIP  optimization  model.  To  set  the  stage  for  the  model 

formulation,  we  first  list  the  eomplete  speeifieations  of  the  predesignation  problem,  as 

addressed  in  the  remainder  of  the  thesis. 

A,  PROBLEM  SPECIFICATIONS:  DEFINITIONS  AND  ASSUMPTIONS 

1,  Near-Simultaneous  Tasks 

•  Eaeh  task  has  an  earliest  time  and  a  latest  time  to  launeh  (ETL,  EXE). 

•  An  average  time  to  launch  (ATE)  (ealeulated  as  the  midpoint  between  ETL  and  ETL) 
ean  be  used  as  hypothetieal  time  to  determine  near-simultaneous  (also  ealled 
conflicting)  tasks. 

•  If  the  differenee  in  the  ATE  for  two  given  tasks  is  less  than  the  given  conflicting  time 
window  (typieally  45  minutes),  both  tasks  must  belong  to  a  eommon  near- 
simultaneous  task  set. 

•  A  task  may  be  in  more  than  one  near-simultaneous  task  set. 

2,  Firing  Units 

•  There  are  two  types  of firing  units:  ships  and  submarines. 

•  Some  firing  units  incur  a  so-called  employment  penalty.  (Related  objeetive:  Goal  #2.) 

•  Eaeh  firing  unit  is  defined  as  either  expend  or  non-expend.  It  is  preferable  to  use 
weapons  from  expend  firing  units.  (Related  objeetive;  Goal  #3.) 


FIRING  UNITS 

1 

f 

▼ 

SHIPS 

1 

SUBMARINES 

Half-Modules 

1  ^  1 

(Weapons) 

Torpedo  Room  Canister  Capsules 

(Weapons)  (Weapons) 

Eigure  10.  Eiring  Units. 

There  are  two  types  of  firing  units:  Ships  and  submarines.  Ships  stores  weapons  in  eells 
in  half-modules,  whereas  submarines  have  two  types  of  storage  for  weapons:  torpedo 
room  and  eanister  eapsules. 
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3,  Ships 

Each  ship  has  a  given  number  of  half-modules,  depending  on  ship  type  as  shown  in 
Figure  2.  Eaeh  half-module  may  eontain  up  to  four  TLAMs,  and  these  ean  be  of 
different  type  (see  weapon  types  below). 

The  number  of  TLAMs  fired  from  a  ship  is  restrieted  by  a  power  eonstraint,  whieh 
allows  only  one  TEAM  to  be  aligned  per  half-module  within  a  eonflieting  time 
window. 

4,  Submarines 

Each  submarine  ean  fire  TLAMs  from  its  torpedo  tubes  and  from  its  eanister  eapsules 
(if  any).  The  number  of  tubes  and  eapsules  depends  on  the  submarine  type. 

Torpedo  tubes  ean  fire  the  TLAMs  available  in  the  torpedo  room.  The  number  and 
type  of  weapons  in  the  torpedo  room  are  submarine  dependent. 

The  maximum  number  of  weapons  that  ean  be  fired  from  the  torpedo  tubes  per 
eonflieting  time  window  is  preeisely  the  number  of  torpedo  tubes. 

Eaeh  eanister  eapsule  ean  fire  only  one  weapon  (the  one  that  is  already  loaded  in  the 
eapsule,  if  any). 

Canister  eapsules  on  a  submarine  are  not  affeeted  by  eonflieting  time  windows. 

Each  submarine  may  or  may  not  be  designated  for  back-up  task-part  alloeation  by  the 
Tomahawk  Strike  Coordinator  (TSC). 

5,  Targets 

Eaeh  target  eonsists  of  one  or  more  tasks. 

Eaeh  target  is  designated  either  as  restricted  or  unrestricted.  A  restrieted  target  needs 
all  tasks  and  task-parts  to  be  alloeated  to  a  limited  number  of  firing  units,  given 
independently  for  eaeh  restrieted  target.  An  unrestrieted  target  has  no  limitation  in 
this  regard. 

A  target  is  deemed  alloeated  when  all  of  its  tasks  have  been  alloeated.  If  the  target  is 
unrestrieted,  it  is  allowed  (although  not  desirable)  that  some  tasks  in  the  target  are 
alloeated  while  others  are  not.  However,  if  the  target  is  restricted,  then  all  or  none 
tasks  of  the  target  tasks  should  be  alloeated. 

6,  Tasks,  Task-Parts  and  Task-Part  Missiles 

Eaeh  task  belongs  to  a  speeific  target. 

There  are  four  types  of  task-parts:  primary  (PR),  ready-spare  (RS),  back-up  (BU)  and 
ghost  (GH).  (See  List  of  Terms  and  Abbreviations.) 

Eaeh  task  may  have  one  or  more  types  of  task-parts. 

Eaeh  task-part  in  a  task  requires  one  or  more  missiles. 


•  Every  task  must  have  either  a  PR  task-part,  or  a  GH  task-part,  but  not  both.  The 
former  is  ealled  a  regular  task,  and  may  also  include  the  following  types  of  task-parts: 
RS  and/or  BU.  If  the  task  has  a  GH  task-part,  it  is  called  a  ghost  task,  and  it  cannot 
include  any  other  type  of  task-part  (see  Figure  1 1). 

•  For  regular  tasks,  there  cannot  be  more  RS  missiles  required  than  PR  missiles 
required. 

•  A  task  is  considered  allocated  only  when  all  of  its  task-parts  have  been  fully 
allocated,  i.e.,  when  all  the  missiles  required  for  each  of  its  task-parts  have  been 
allocated.  Otherwise,  the  task  is  deemed  unallocated.  (Partial  allocation  of  a  task  is 
not  allowed.)  (Related  objective:  Goal  #1.) 

•  PR  and  RS  task-parts  must  be  allocated  to  the  same  firing  unit  in  pairs:  PRl-RSl 
first;  then  PR2-RS2,  and  so  forth.  The  nth  primary  task-part  (PRn)  can  be  allocated 
alone  if  RSn  is  not  required. 

•  There  is  no  limit  in  the  number  of  missiles  required  for  a  BFl  task-part  in  a  regular 
task,  nor  it  is  required  that  a  RS  task-part  exists  for  the  task  in  order  to  require  one  or 
several  BU  task-parts. 

•  BU  task-parts  (for  a  given  task)  must  be  allocated  to  firing  units  excluding  all  of  the 
firing  units  used  for  the  PR  task-parts  (for  the  same  task). 

•  Each  task-part  of  each  task  is  classified  as  separable  or  non-separable.  All  the  missile 
requirements  of  a  non-separable  task-part  must  be  allocated  to  the  same  firing  unit. 
(Separable  and  Non-separable  are  sometimes  referred  to  as  split  and  non-split, 
respectively) 


19 


7.  Launch  Areas 

•  Each  firing  unit  is  assumed  to  be  loeated  in  a  unique  launch  area,  and  eannot  be 
moved  to  a  different  area. 

•  Eaeh  launeh  area  has  a  list  of  targets  that  can  be  attaeked  from  the  area.  A  firing  unit 
in  the  launeh  area  of  a  target  is  eonsidered  geo-feasible. 

•  Primary  task-parts  alloeated  to  firing  units  in  a  partieular  area  ean  be  eharaeterized  as 
spread  or  not  spread.  This  refers  to  whether  or  not  spreading  the  primary  task-parts 
aeross  the  different  firing  units  in  the  area  is  desired.  (Related  objeetive:  Goal  #5.) 

•  Baek-up  task-parts  alloeated  to  firing  units  in  a  partieular  area  ean  also  be 
eharaeterized  as  spread  or  not  spread.  (Related  objeetive:  Goal  #6.) 
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8,  Weapons 

•  Each  missile  required  by  any  task-part  in  a  given  task  must  be  allocated  a  weapon 
type  from  the  so-ealled  (missile  to  mission  matching)  list.  This  list  is  task- 
dependent. 

•  Each  available  weapon  in  the  M  list  for  the  task  has  an  assoeiated  priority  to  be  used 
that  makes  it  more  or  less  desirable  to  use  compared  to  the  other  weapons  on  the  list. 
(Related  objective:  Goal  #7.) 

•  Each  weapon  type  has  a  value  that  is  used  to  compute  the  residual  salvo  value  after 
allocation.  Residual  salvo  computation  must  take  into  account  half-module  and 
torpedo  tube  power-constraints.  Thus,  only  the  most  valuable  weapon  from  each 
half-module,  and  a  limited  number  of  weapons  from  the  torpedo  room  (as  many  as 
torpedo  tubes)  must  eount  for  the  residual  salvo  computation.  (Related  objective: 
Goal  #8.) 

9,  Tomahawk  Strike  Coordinator  Allocation 

The  TSC  may  specify  mandatory  allocation  of  selected  task-part  missile 
requirements  by  providing,  for  each  of  these  requirements,  the  following  information: 

•  A  ship,  a  half-module,  a  weapon  type  (for  alloeation  to  ships),  or 

•  A  submarine  and  either  a  weapon  type  in  the  torpedo  room  or  a  weapon  type  in 
eanister  eapsules  (for  allocation  to  torpedo  tubes  or  canister  capsules  on  submarines, 
respeetively). 

10,  Goal  Specifications 

a.  Goal  #1 

“Meet  All  Assigned  Tasking”  -  To  date,  this  goal  has  been  addressed  by 

minimizing  the  number  of  “unassigned  task-part  missile  requirements.”  This 
interpretation  is  assumed  in  our  formulation. 

In  some  cases,  however,  it  might  be  desirable  to  address  Goal  #1  by 
minimizing  some  of  the  following: 

•  A  weighed  average  of  unmet  tasks-parts.  Eor  example,  eonsider  the  following  task- 
part  missile  requirements: 

•  Task  1 :  PRl ,  BUI ,  BU2,  BUS 

•  Task2:  PR1,PR2,  PR3 


Sinee  a  task  eannot  be  alloeated  unless  all  of  its  task-parts  are  alloeated. 


eurrently  it  is  deemed  more  important  to  alloeate  the  four  missiles  required  for  Task  1 
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(even  if  three  of  them  are  BUs)  than  the  three  missiles  required  for  Task  2,  all  of  whieh 
are  PRs.  This  ean  be  remedied  by  weighing  the  task-parts. 

•  The  number  of  unmet  tasks,  instead  of  the  number  of  unmet  task-parts.  These  two 
outputs  might  be  signifieantly  different,  depending  on  how  evenly  the  missile 
requirements  are  spread  among  the  tasks. 

•  The  number  of  targets  that  eannot  be  fully  alloeated.  (Non-restricted  targets  allow 
some  tasks  to  be  alloeated  even  if  other  tasks  in  the  target  eannot  be  alloeated.) 

b.  Goal  #2 

“Minimize  the  Use  of  Firing  Units  Currently  Engaged  in  Other  Tasking”  - 
Currently,  this  goal  has  been  aeeomplished  with  the  use  of  employment  penalties. 

c.  Goal  #3 

“Maximize  Missiles  Alloeated  from  Designated  Expend  Eiring  Units”  - 
Remark:  Only  alloeated  missiles  for  PR  task-parts  are  deemed  used.  That  is,  even  if 
other  RS,  BU  and  GH  missiles  are  alloeated,  these  do  not  beeome  part  of  the  preliminary 
eount  of  fired  missiles. 

d.  Goal  #4 

“Eevel  the  Number  of  Missiles  Remaining  on  Non-Expend  Eiring  Units”  - 
Remarks:  Unalloeated  missiles  along  with  those  missiles  alloeated  to  RS,  BU  and  GH 
task-parts  are  deemed  remaining.  The  level  definition  is  based  on  the  total  differenee 
between  the  arithmetie  average  of  the  remaining  missiles  aeross  all  non-expend  firing 
units  and  the  remaining  missiles  on  eaeh  non-expend  firing  unit  (even  if  it  is  a  firing  unit 
with  positive  employment  penalty). 

e.  Goal  #5 

“Spread  Primary  Task-Part  Missiles  Among  As  Many  Eiring  Units  As 
Possible,  in  Eauneh  Areas  So  Designated.” 

f.  Goal  #6 

“Spread  baek-up  task-part  missiles  among  as  many  firing  units  as 
possible,  in  launeh  areas  so  designated.” 

g.  Goal  #7 

“Prom  among  all  feasible  missile  types  for  eaeh  task,  ehoose  the  least 
eapable  missile,  so  as  to  preserve  more  eapable  missiles  for  the  future.” 
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h.  Goal  #8 

“Maximize  residual  salvo  capability”  -  Remark:  Only  allocated  PR  task- 
parts  are  used  to  compute  the  post-allocation  salvo  value.  To  do  this  for  ships,  we  use  the 
most  valuable  weapon  in  each  half-module.  For  submarines,  we  use  the  most  valuable 
weapons  in  the  torpedo  room  that  can  be  fired  simultaneously  from  the  torpedo  tubes,  and 
all  the  weapons  in  the  canister  capsules. 

B,  TLAM  ALLOCATION  MODEL 

In  this  section,  we  present  the  problem  formulation.  In  order  to  simplify  the 


presentation,  we  have  divided  the  formulation  into  several  parts  driven  by  the  goal  and/or 

the  feature  in  the  specifications:  Initially,  we  present  a  basic  formulation  for  the  first  goal 

with  limited  problem  features.  Then,  we  add  other  goals  and  specifications  such  as 

submarines. 

restricted  targets,  TSC  specifications,  etc.  The  full  Allocation  Model  is 

shown  in  Appendix  A. 

1. 

Basic  Model  with  Goal  #1 

In  this  model. 

we  only  consider  ships  as  firing  units,  and  goal  #1  in  the  objective 

function. 

•  Indices 

set  of  ships;  s  eS 

H, 

set  of  half-modules;  h  eH  =  {1,2,3,...} 

G, 

set  of  targets;  g  e  G 

T, 

set  of  tasks;  t  eT 

P, 

set  of  task-parts;  p^P  =  [PR,RS,BU  ,GH^ 

M, 

set  of  task-part  missile  requirements;  m,m'  eM  =  {1,2,3,...} 

w, 

set  of  weapon  types;  w^W 

I, 

set  of  subsets  of  near-simultaneous  tasks;  i  e  / 

A, 

set  of  launch  areas;  A 

K, 

set  of  goals;  A:  e  A  =  {l,2,...,8}  .  (Remark:  We  introduce  the  set  of 

all  goals  although  we  only  formulate  goal  #1  here.  The  rest  of  the 
goals  will  be  formulated  later.) 
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•  Sets 


Tjf ,  subset  of  tasks  for  target  g;  cz  T 

,  subset  of  near-simultaneous  tasks  in  group  i\  T‘  a  T 

Gf,  subset  of  targets  whieh  are  geo-feasible  from  area  a; 

Wj ,  subset  of  weapons  that  ean  be  used  for  task  t\  Wj  c  W 


•  Data 


^shw^ 

H 


tp  ■ 


number  of  weapons  of  type  w  in  half-module  h  on  ship  s 

number  of  half-modules  on  ship  s 

number  of  missiles  required  for  task-part  p  of  task  t 

’  area  for  ship  s 

•  Derived  Data  and  Sets 

g^,  target  for  task  t  (i.e.,  t 

Sf ,  subset  of  geo-feasible  ships  for  task  t;  Sf  =  is  g^  eG 


•  Deeision  Variables 

shwtpm’  equals  1  if  a  weapon  of  type  w  in  half-module  h  of  ship  s  is 
used  to  alloeate  missile  m  of  task-part  p  in  task  t,  0  otherwise 

yS 

stpm  ’  equals  1  if  missile  m  of  task-part  p  in  task  t  is  alloeated  to  ship 
5,  0  otherwise 

^tpm’  equals  1  if  task  t  is  unmet,  0  otherwise 
Z 

objeetive  funetion  value  for  goal  k^K.  (Remark:  We 
introduee  deeision  variables  for  all  objeetive  funetions  here 
although  we  only  use  goal  #1  in  this  model). 


•  Formulation 

Goal  k=  \  :  Minimize  the  total  number  of  unmet  task-parts: 
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Z;  =  minZi,  where  =X  Z  Z 


t  p  m=\ 


Subjeet  to: 

-  Ensure  weapon  availability  for  tasks  alloeated  to  half-modules: 


U.  sss  ^ shwtpm  ^shw  ’ 


ys,h,w\<h<m^  ,n,>0 


sgS;,  P  m=\ 
w=Wl 


Establish  the  alloeation  at  the  firing  unit  level: 


1.2.  S  S  ^shwipm-^stpm’  ys,t,p,m\s  ,l<m<n^ 


h=\  WGfffnL>0 


1  ^  y  7^  <1 

t.i.  ^  ^stpm  — 


1  <  m<  n 


ses; 


Ensure  all  (or  none)  of  the  task-parts  in  a  task  are  alloeated: 

\/f  rt  m  ll  <  m  < 

14  y  ,  =  y  7 

/  V  stp  m  /  j  s 


tS 

'■  stpm 


\/t,p,ml<m<n^  ,n^,  >  0, 


.es/  ses;  for  p'^{'PR","GH"},m'  =  "r 

-  Prevent  near-simultaneous  tasks  in  the  same  half-module: 


1.5.  Z  Z 


ys,h,i 


l<h<m 


H 


Alloeate  eaeh  task-part  missile  requirement  or  take  a  penalty: 


i-e-  Z  C+Vp,.  =1. 

s&Sj 


yt,p,m  \  <m<n^ 


1  Y  ys  I  y5  <1 

^  •  ^stpm  stp'ni  —  ^ 


Ensure  primary  and  baek-up  task-part  relations: 

ys,t,m,m'\^s  ey  ,1  <  m  <  1  ^  m  <n^,, 

for  p  =  "PR",p  ="BU" 

Ensure  primary  and  ready-spare  task-part  relations: 


1  8 

’  *  ^  Stpm  ^stp'm^ 


for  p  =  "PR",p="RS" 
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-  Binary  variables: 


\fs,h,w,t,p,ms  ,\<h<mf  ,weJ^  ,t  eT,r(f^>0,peP,\<m<n^ 


-  Non-negative  variables  and  bounds: 

1.10.  0<7j^^<l,  '\fs,t,p,mseSf,teT,peP,l<m<n^ 

1.11.  0<17,^„<1,  yt,p,mteT,peP,l<m<n^ 


2.  Goal  #2 

•  Data 

p^,  employment  penalty  for  ship  5 

•  Deeision  Variables 

,  equals  1  if  ship  5  is  engaged  in  any  tasking  in  target  g,  0  otherwise 
,  equals  1  if  ship  5  is  used,  0  otherwise 

•  Formulation 

Goal  k  =  2:  Minimize  the  use  of  penalty  firing  units: 

Z’  =  minZ2,  where  =  '^p^Z^ 

s 

Additional  eonstraints: 


-  Caleulate  firing  unit  engagement  in  targets  and  in  the  overall  alloeation: 


2  1  Y  <ZZ 

^stpm 


y s,t, p,m,g  s  e  Sf  ,\<  m<  n 


M 

ip 


2.2.  Z4<Zf,  V5,g 

(Remark:  Engagement  is  ealeulated  at  the  level  of  targets  for  its  use  in  other 
speeifieations  later  in  the  formulation.  For  the  purpose  of  ealeulating  firing  unit 
engagement  in  the  alloeation  only,  (2.1)  and  (2.2)  ean  be  eombined  into: 
Ysl„  <Z^,ys,t,p,m\seSf,\<m< n" ) 

-  Binary  variables: 

2.3.  ZZjejO,!},  V5,g|5  e  5',g  e  G 
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2.4.  Zfe{0,l}, 


\/s  eS 


3.  Goal  #3 

•  Sets 

subset  of  expend  ships;  5*^^^  c  S 

•  Decision  Variables 

,  total  number  of  weapons  allocated  to  primary  task-parts  from 
expend  ships 

•  Formulation 

Goal  k  =  3:  Maximize  the  number  of  weapons  allocated  from  designated  expend 
firing  units: 

Z3  =  maxZj,  where 

Additional  constraints: 

-  Use  of  expend  firing  units  for  primary  task-parts: 

3.1.  z  Z  Z  C.  for  p  =  "PR" 

t^SESj'  "'=1 

-  Non-negative  variables: 

3.2.  E^>0 

4.  Goal  #4 

•  Data 

nNE.p^ 

•  Decision  Variables 

RN^, 

RAVG, 

RD^, 


number  of  non-expend  firing  units:  =  p  \  5*^^^  | 

number  of  missiles  on  non-expend  ship  s  that  are  not  allocated 
to  a  primary  task-part 

average  of  the  residual  number  of  missiles  across  non-expend 
firing  units  after  allocation 

difference  between  the  residual  number  of  missiles  on  ship  s 
and  RA  VG 
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•  Formulation 

Goal  k  =  A:  Minimize  the  deviation  from  average  of  residual  number  of  missiles 
on  non-expend  firing  units; 

Z*  =  minZ4,  where  Z^=  ^ 

Additional  eonstraints: 

-  Compute  absolute  deviation  from  average  number  of  residual  weapons  onboard 
each  firing  unit: 

4-1-  Z  ZoL-Z  Z  C  lor  p  =  "PR" 

h=l  w  tlsi^Sj  «i=l 


4.2.  ^  RN^  =n^^^‘’RAVG 


4.3. 

RN^  -RAVG>-RD^, 

VseS\S^^^ 

4.4. 

RN^  - RAVG  <  RD^ , 

VseS\S^^^ 

-  Non-negative  variables: 

4.5. 

RN^  >  0, 

Vs  e  S 

4.6. 

IV 

o 

Vs  e  S 

4.7. 

RAVG>0 

5.  Goal  #5 

•  Sets 

Z™,  subset  of  areas  whose  primary  task-parts  are  desired  to  be  assigned 
to  as  many  firing  units  as  possible;  A™  a  A 

•  Decision  Variables 

PR^ ,  equals  1  if  ship  5  is  assigned  to  allocated  any  primary  task-part,  0 
otherwise 

•  Formulation 

Goal  k  =  5:  Maximize  the  number  of  firing  units  that  have  been  assigned  to 
allocate  primary  task-parts  in  launch  areas  so  designated: 
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6.  Goal  #6 

•  Sets 

,  subset  of  areas  whose  back-up  task-parts  are  desired  to  be  assigned 
to  as  many  tiring  units  as  possible;  a  A 

•  Decision  Variables 

,  equals  1  if  ship  5  is  assigned  to  any  back-up  task-part,  0  otherwise 

•  Formulation 

Goal  k  =  6:  Maximize  the  number  of  firing  units  that  have  been  assigned  to 
allocate  back-up  task-parts  in  launch  areas  so  designated: 


Z;  =  maxZg,  where  Z,  =  Y,  Z 

aEBU'‘  \ys\a^=a  ^ 

Additional  constraints: 

-  Keep  track  of  areas  that  must  spread  primary  task-part  missiles: 

(>■1.  S  E  t0Tp  =  "BU 

>n=l 

-  Binary  variables: 

6.2.  5f/f  e{0,l},  VseS 
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7.  Goal  #7 

•  Data 

-5  ^ 

wip ,  weight  of  task-part  p  to  match  the  M  list  position  of  for  task-part  p 
m]^,  position  of  weapon  w  e  for  task  t 

•  Formulation 

■y 

Goal  k  =  1:  Minimize  the  missile  to  mission  matching  (M  )  list  position  of  all 
selected  weapons: 

Z;  =  min  Z„  where  ^  ZZZ 

a  h=\  w|Bfj„>0  IseS/’,  P 
\weW/' 


8.  Goal  #8 

•  Data 

v^,  relative  value  of  weapon  w  for  residual  salvo  computation 

•  Decision  Variables 

D^hw’  equals  1  if  one  or  more  weapons  of  type  w  remain  on  ship  5  in  half¬ 
module  h  after  firing  all  primary  task-part  missiles,  0  otherwise 

•  Formulation 

Goal  k=  S:  Maximize  residual  salvo  capability: 

mf 

Z;  =  maxZ8,  where  Z,  =  ^  ^  ^ 


Additional  constraints: 

-  Calculate  residual  salvo  size: 

8 •  1  •  -  Z  Z  ^sh^tpm  ^  aL .  ys,h,w 

lieSf,  m=l 

8-2.  Z  AL  ys,h\\<  h  <  m^ 

w|nV>0 

-  Non-negative  variables  and  bounds: 

8.3.  0<Df^„<l,  \/s,h,w\s  eS,\<h<mf  ,weW,'^  ,n^f,„>0 


\l<h<mi 


for  p  =  ''PR'' 
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9. 

•  Sets 


Restricted  Targets 


G*,  subset  of  targets  whieh  are  designated  as  restricted',  (^G 


•  Data 


•  Formulation 


maximum  number  of  firing  units  to  whieh  a  restrieted  target  g  can 
be  alloeated 


M 


Additional  eonstraints; 

-  Caleulate  firing  unit  engagement  in  the  target: 

Same  as  eonstraint  (2.1). 

stpm  ^  t,p,m,g\s&Sj  ,\<m< 

-  Limit  firing  units  for  restrieted  targets: 

9.2.  ^ 

-  Alloeate  all  (or  none)  of  the  tasks  in  restrieted  targets: 


0  1  <77^ 

^stDm  —  ’ 


93  u  =u , 

^  tpm  ^  t'pm^ 


'dt,t',g 


n"  >0,for  p^{'PR'\"GH"],m  =  "\" 
t,t'  eT^,t^t',geG'^ 


10,  Non-separable  Tasking 

•  Sets 

,  subset  of  task-parts  in  task  t  that  are  not  separable;  a  P 

•  Deeision  Variables 

NSlp  ,  equals  1  if  all  missile  requirements  for  non-separable  task-part  p  of 
task  t  are  assigned  to  ship  5,  0  otherwise 

•  Formulation 

Additional  eonstraints: 

-  Ensure  the  alloeation  of  all  missiles  for  a  non-separable  task-part  to  the  same 
firing  unit: 
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11,  TSC  Specifications 

•  Sets 

TSC^ ,  subset  of  six-tuples  {s,h,w,t,p,m)  sueh  that  the  TSC  speeifies  that 
task-part  missile  {t,p,m)  must  be  alloeated  to  weapon  type  w  in 
half-module  h  on  ship  s. 

•  Formulation 

Additional  eonstraints: 


-  TSC  speeifioations; 


11.1. 


^  shwtpm 


=  1, 


V {s,  h,w,  t,  p,  ni)  eTSCf 


12,  Submarines 

Most  sets  and  data  for  ships  have  an  equivalent  notation  for  submarines.  Model 
eonstraints  ean  be  extended  in  order  to  ineorporate  submarines.  Some  of  these  extensions 
are  immediate.  Others,  however,  require  more  elaboration  due  to  differenees  in 
speeifioations  for  ships  and  submarines.  Here  is  an  example; 

•  Indioes 

B,  set  of  submarines;  b^B 

•  Sets 

Bf ,  subset  of  geo-feasible  submarines  for  task  t;  Bf  cz  B 

•  Data 

,  number  of  weapons  of  type  w  in  the  torpedo  room  of  submarine  b 
number  of  torpedo  tubes  on  submarine  b 

•  Deoision  Variables 

X^l^pm,  equals  1  if  a  weapon  of  type  w  in  the  torpedo  room  of 

submarine  b  is  used  to  alloeate  missile  m  of  task-part  p  in  task 
t,  0  otherwise 
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•  Formulation 


Additional  constraints: 

-  Avoid  near-simultaneous  task  alloeation  to  torpedo  tubes  (similar  to  (1.5)  for 


ships): 


'V 


12.1.  X  X 


teT 


\/b,i 


6gS;,  P  m=l 


Binary  variables 


12.2.  A-".€{0,1}, 


Vb,  w,t,p,m  b  e  B,w€  WJ  ,  ufj  >Q,t  eT  ,p  e  P,\<m<nl 


Aeeordingly,  the  formulation  for  goals  and  other  eonstraints  must  be  extended  in 
order  to  ineorporate  submarines.  The  eomplete  formulation  is  shown  in  the  Appendix. 

13.  Levels  of  Achievement  for  Goals 

When  goal  A:  >  1  is  being  optimized,  it  is  possible  to  speeify  a  target  level  for  all 
higher  priority  goals,  k'  =  \,...,k -\ . 

•  Data 

fj,  pereentage  of  deviation  (from  the  best  possible)  allowed  for  goal  k 
when  a  less  priority  goal  is  eonsidered  in  the  objeetive  funetion 
direetion  of  goal  k:  ^^.=1  if  maximization,  ^<-=-1  if 
minimization 

•  Formulation 

Additional  eonstraints: 

-  Foree  levels  of  aehievement  for  goals: 

13.1.  z,,<zl,{\+s,),  \fk'\k'<k,  r,  =  -i 

13.2.  Z,,>Zl,{\-s,),  \fk'\k'<k,  r,.=\ 

14,  Combined  Objective  Function 

As  an  alternative  to  hierarehieal  optimization,  we  may  optimize  a  eombined 
(weighed)  objeetive  funetion  of  all  goals. 
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•  Data 


a>i^,  weight  given  to  goal  k  when  a  eombined-goal  objeetive  is 
considered,  where  <»,  »  Oj  ^  ^  • 

•  Decision  Variables 

Z^,  combined-goal  objective 

•  Formulation 

Additional  constraints; 

-  Combined  goal;  Minimize  a  weighed  average  of  the  goals  k  e  K  : 

14.1.  Z*  =  minZo,  where 

k 

15.  Multi-Objective  Formulation 

Appendix  A  contains  one  formulation  (A.1)-(A.73)  that  is  the  union  of  all  the 
previous  formulations.  We  consider  two  formulations  for  our  multi-objective  problem; 

•  Hierarchical  formulation  (H);  For  each  k  e  K ,  solve  the  following  sequence  of 
models; 

(//^);  Z^  =  min  (or  max)  Z^ 
s.t.  (A.1)-(A.73) 

Remark;  To  be  more  precise,  we  write  all  these  models  as  the  following 
minimization  problems; 

(Hk)-  Zl=-r,  min- 

s.t.  (A.1)-(A.73) 

•  Combined-goal  (sometimes  referred  to  as  weighed-sum)  formulation  (C); 

(C);Z:  =  minZ„ 

s.t.  (A.1)-(A.45) 

(A.48)  -  (A.73) 

C.  PRELIMINARY  VALIDATION  OF  THE  MODEL 

In  this  section,  we  test  our  formulation  with  a  small  case  from  “NFS  Report  on 

TLAM  Heuristic”  [2003]  in  which  the  heuristic  algorithm  fails  to  obtain  the  optimal 

solution.  As  we  will  see,  our  model  provides  the  optimal  solution  by  using  both  the 
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Hierarchical  and  Combined-goal  formulations.  Figures  12  and  13  show  the  task 
specifications  and  firing  units  for  the  case,  respectively. 


TASKS  SPECIFICATIONS 


Target 

Task 

Conflict  Set 

#ofPR 

#ofRS 

#ofBU 

#ofGH 

T1 

T1-1 

1 

2 

1 

1 

0 

12 

T2-1 

1 

1 

0 

0 

0 

T2-2 

2 

1 

0 

0 

0 

T2-3 

2 

1 

0 

0 

0 

IIBrilllBJ 

Figure  12.  Task  Specifications  for  Validation  Case. 

We  have  two  targets  and  four  tasks.  Tasks  “Tl-1”  and  “T2-1”  are  near-simultaneous  and 
the  same  occurs  with  tasks  “T2-2”  and  “T2-3.”  The  list  indicates  the  type  of  weapons 
that  can  be  used  for  each  task. 


FIRING  UNITS  (Ships) 

FU1 

H1 

II-C 

II-C 

H2 

II-C 

II-C 

I _ BJ2 _ I 

H1 

H2 

H3 

lll-C 

lll-C 

II-C 

II-C  I 

Figure  13.  Firing  Units  for  Validation  Case. 

There  are  two  ships:  The  first  one  has  two  half-modules  with  III-C  and  II-C  TLAMs, 
while  the  second  one  has  three  half-modules. 


The  solution  provided  by  the  heuristic  algorithm  (as  reported  in  “NFS  report  on 
TLAM  heuristic”  [2003])  is  displayed  in  Figure  14.  It  leaves  task  “Tl-1”  unmet. 


I  NON-OPTIMAL  ALLOCATION  | 

Target 

Task 

Task-Part 

Firing  Unit 

1  Half-Module  | 

Weapon  | 

T1 

T1-1 

PR1 

UNMET 

PR2 

UNMET 

RSI 

UNMET 

BU1 

UNMET 

T2 

T2-1 

PR1 

FU2 

H3 

T2-2 

PR1 

FU2 

H3 

T2-3 

PR1 

FU2 

H1 

lll-C  I 

Figure  14.  Non-Optimal  Solution  of  Validation  Case. 

The  heuristic  solution  cannot  allocate  task  “Tl-1.” 
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Two  solutions  to  this  problem  are  displayed  in  Figure  15  and  Figure  16,  aehieved 
by  the  Hierarehieal  formulation  and  the  Combined-goal  formulation,  respeetively.  As  we 
will  see,  the  first  one  is  optimal,  whereas  the  seeond  one  is  near-optimal. 


OPTIMAL  SOLUTION 

Target 

Task 

Task-Part 

Firing  Unit 

Half-Module 

Weapon  | 

T1 

T1-1 

PR1 

FU2 

H1 

PR2 

FU2 

H3 

RS1 

FU2 

H2 

nsB^m 

BU1 

FU1 

H2 

T2 

T2-1 

PR1 

FU1 

H1 

ll-C  I 

T2-2 

PR1 

FU2 

H1 

isB^m 

T2-3 

PR1 

FU2 

H3 

ll-C 

Figure  15.  Optimal  Solution  of  Validation  Case  Obtained  with  the  Hierarehieal 

Formulation. 


OPTIMAL  SOLUTION 

Target 

Task 

Task-Part 

Firing  Unit 

Half-Module 

Weapon  | 

T1 

T1-1 

PR1 

FU2 

H1 

PR2 

FU2 

H2 

isa^m 

RS1 

FU2 

H3 

BU1 

FU1 

H1 

isa^m 

T2 

T2-1 

PR1 

FU1 

H2 

T2-2 

PR1 

FU2 

H2 

isa^m 

T2-3 

PR1 

FU2 

H3 

ll-C 

Figure  16.  Near-Optimal  Solution  of  Validation  Case  Obtained  with  the  Combined- 

Goal  Formulation. 

We  have  used  the  following  weights:  wi=  10^,  W2=  10^,  ...,  W8=  1. 


Complete  lists  of  goal  aehievements  for  both  formulations  are  shown  in  Figure 
17.  (Remark:  The  Hierarehieal  formulation  yields  a  better  solution  than  the  Combined- 
goal  formulation  for  Goal  #8  only) 
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Goal 

Priority 

Goal  Name 

Direction 

Hierarchical-goal 

Formulation 

Combined-goal 

Formulation 

1 

Unmet  Task-parts 

Min 

0.0 

0.0 

2 

Employment  Penalty 

Min 

0.0 

0.0 

3 

Expend  Firing  Units 

Max 

0.0 

0.0 

4 

Level  Firing  Units 

Min 

1.0 

1.0 

5 

Primary  Spread 

Max 

0.0 

0.0 

6 

Back-up  Spread 

Max 

0.0 

0.0 

7 

List  Position 

Min 

25.0 

25.0 

8 

Residual  Salvo 

Max 

10.0 

0.0 

Figure  17.  Goal  Achievements  for  Validation  Case. 

Both  Hierarchical  and  Combined-goal  formulations  find  satisfying  solutions,  but  the 
solution  obtained  by  the  Hierarchical  formulation  is  superior  to  that  of  the  Combined- 
goal  formulation  for  the  “Salvo  Capability”  goal. 

Further  validation  of  our  formulations  is  extensively  analyzed  in  Section  Ill. 

D,  SOLVING  TLAM  ALLOCATION  MODEL  USING  FIX-AND-RELAX 

In  this  section,  we  describe  our  use  of  a  general-purpose  methodology,  Fix-and- 
Relax  (F&R  )  [Dillenberger  et  ah,  1994],  to  tackle  some  of  the  large-scale  cases  that  we 
present  in  Section  III.  (Similar  techniques  were  used  by  Brown  et  al.  [1987]  under  the 
name  of  “Cascading.”)  Branch-and-Bound  (B&B)  [Wolsey,  1998]  can  become 
inefficient  to  solve  MIPs  when  the  numbers  of  integer  variables  is  large.  In  order  to 
reduce  the  computational  burden  of  solving  our  model  with  B&B,  we  apply  F&R  ,  which 
involves  solving  a  number  of  sub-problems,  each  one  with  smaller  complexity  than  the 
original  problem.  In  our  model,  we  do  not  implement  F&R  for  all  of  the  binary 
variables,  only  those  from  which  the  other  binary  variables  are  derived  e  {0,l}, 

and 

For  example,  using  the  combined  goal  formulation,  problem  (C)  can  be  rewritten 
as  follows: 
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(C):min/(x,j,M) 
[xe{0,l}^ 


s.t. 


M  >  0 

x,y,u  e  XYU 


where  x  is  an  A^-dimensional  veetor  comprised  of  (see  Appendix  A): 

^  ,  for  s  eS,\<h<m^  ,wewj'  ,t  eT,n^j  >0,peP,l<m<n 

shwtpm  s  t  shw 

^bwtpm’  ^  ^  >Q,t &T,p &P,\<m<  n",  and 

^bwtpm^  b^B,w^W^,nZ>0,t^T,p^P,\<m<nf^, 


M 

tp  ’ 


while  y  is  the  M-dimensional  vector  of  the  remaining  binary  variables,  which  appear  in 
constraints  (A.48)-(A.60)  in  the  original  model.  In  this  model,  u  is  the  vector  of  all 
continuous  variables,  which  appear  in  constraints  (A.61)-(A.73).  f{x,y,u)  is  our  linear 
objective  function.  The  constraints  represented  by  x,y,u  ^  XYU  are  (A.1)-(A.45)  and 
(A.48)-(A.73)  for  the  (C)  model  or  (A.1)-(A.73)  for  each  of  the  (77^^)  models. 


Let  Xj,...,x„  be  the  components  of  x  where  N  is  the  total  number  of  binary 
variables  to  which  we  implement  F&R ,  that  is,  Xj  in  our  abbreviated  formulation  is 
representing  a  variable  of  type  or  XZp„  or  .  Let  L  =  {l,....  A}  be  the  set  of 

indices  for  these  variables,  and  be  a  partition  of  the  set  V.  That  is, 

n 

V.  ^  V,Vi  =  and  V.  nV.r  =0,\/iJ' =  ^  i'.  The  cardinality  of 

i=\ 

n 

each  V.  is  denoted  |fj.  |  =  A,. ,  so  N  =  . 

i=i 
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Now,  we  rewrite  our  (C)  model  as: 


(C):  min/(x,j,M) 

e{0,l},y/'el^,z  =l,...,n 

M  >  0 

x,y,u  e  XYU 

F&R  approximates  the  solution  to  this  problem  by  solving  a  sequenee  of  n 
mixed-0-1  sub-problems  (hereafter  ealled  denoted  as  (O)  ,  for  r  =  \,...,n.  (C^) 

ean  be  written  as  follows: 

(C):  min  f{x,y,u) 

{x,y,u)^XYU 

w>0 

Xj=x.,yj  ^VJ  =  \,...,r-\ 
s.t.<  x^-  e{0,l},y/'eF^ 

x^  e[0,l],y7'el^.,z  =  r + 

where  the  values  Xj  for  j  ^V.,i  =  \,...,r-\  in  stage  r>\  are  those  obtained  from  the 

solution  of  problems  (C'),...,  (C^"')  ,  respeetively.  (C^)  ean  be  solved  effieiently,  sinee 
only  a  redueed  subset  of  0-1  integer  variables  must  be  determined  at  eaeh  stage  r  . 

We  seleet  our  variable  partition  by  grouping  tasks.  We  ereate  n  groups  of  tasks, 

n 

I|  ,  sueh  that  \^T  .=  T  .  FJ  is  the  set  of  variables  indiees  assoeiated  with  tasks  in 

/=1 

T, ,  that  is,  with  variables  ,  X^pm  ,  for  t  e  7; . 

We  begin  our  implementation  by  relaxing  the  binary  eonstraints  for  the 
X- variables  exeluding  the  ones  assoeiated  with  tasks  in  7j,  whieh  are  foreed  to  be 

binary.  After  the  initial  solution  obtained  from  (C^) ,  these  binary  variables  are  fixed  in 
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(C^) .  Then,  we  foree  a  new  set  of  variables  (assoeiated  with  tasks  in  T^)  to  be  integer, 
and  repeat  the  proeess  until  we  solve  the  last  stage  model  (C”) . 


Model  (C^)  for  r  =  ean  be  stated  as: 

(C^) :  min  / {x,y,u) 

'(A.l)  to  (A.45)  and  (A.48)  to  (A.73) 

X^,  =X^,  ,  for  s  e  S  ,l  <  h  <  ,w  e  ,n^j  >  0,  p  e  P,l  <  m  <  ,t  e  T.,  i<r 

shwtpm  snwtpm  s  t  snw  tp  ' 


S.t.i 


^BT 

bwti 

■  BC 


BT 

bwti 

-BC 


^bwtpm=^bwtpm’  beB,weW,\n::  >0,peP,l<m<nf^,ter,  i<r 


^bwtpm  ^^bwtpm’  beB,weW/,n:^  >0,peP,l<m<nr^,teT,,  i<r 


shwtpm 

BT  g 
bwtpm 

bwtpm 

s 

shwtpm 

BT  g 
bwtpm 

BC  ^ 
bwtpm 


m^ 

s 

V 

o 

G 

VI 

m  < 

M 

tp 

^BT 
’  ^bw 

V 

o 

m 

P,l  <  m  < 

6  T 

.BC 
’  ^bw 

V 

o 

m 

P,\<  m  <  nfp 

6  T 

m^ 

,w 

>  0,p 

6 

p,l  < 

m  < 

M 

n 

s 

t 

shw 

tp 

■>  '^bw 

V 

o 

m 

S 

VI 

af 

eT, 

i  >  r 

BC 

V 

o 

m 

P,\  <  m 

i  >  r 

Let  Z*(C)  represent  the  optimal  objeetive  funetion  value  for  model  (C) .  Sinee 
(C')  is  a  relaxation  of  (C) ,  the  optimal  objeetive  funetion  value  of  this  model,  Z*(C') 
is  a  lower  bound  on  Z*(C)  .  Similarly,  sinee  (C”)  is  a  restrietion  of  (C),  its  objeetive 
funetion  value,  Z*(C")  is  an  upper  bound  on  Z*(C),  and  provides  us  with  a  feasible 
solution  to  problem  (C)  . 

Basic  F«&:R  Algorithm  for  Model  fC) 

Input:  Partition  Vy,...,V^  ,  where  eaeh  fj  eontains  exaetly  all  the  x -variables 
assoeiated  with  period  tasks  in  set  T. : 
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r,  DT" 

r=  (b.w.I.p.m)  for  and 


,  for  fe?;},  Vi  =  l,...,n 

bwtpm  ‘ 


Step  1:  Set  i  =  1  and  solve  (C‘) 


If  (C‘)  is  infeasible,  STOP:  “Problem  (C)  is  infeasible.”  Otherwise,  set 


^(C)=Z‘(C'). 


Step  2:  If  /  =  n ,  set  Z*(C)=  Z*(C‘)  and  STOP:  “Problem  (C)  is  feasible.” 
Otherwise,  inerease  i  by  1 . 

Step  3:  Solve (C‘)  (see  remark  below).  If  (C‘)  is  feasible,  STOP:  “Problem  (C) 
status  is  unknown.” 

Otherwise,  go  baek  to  Step  2. 

Output:  (C‘)  status  “Infeasible,”  “Feasible”  or  “Unknown.”  If  status  is 

“Feasible,”  ^(C)  and  Z*(C)  are  a  lower  and  an  upper  bound, 
respeetively,  on  the  optimal  solution  to  (C)  . 

Remark:  If  all  variables  in  V-  are  integer  after  solving  (C‘''),  then  this  solution  is 

already  optimal  to  (C‘) 

This  basie  F&R  algorithm  eould  fail  to  provide  an  answer  to  our  problem  if  (C^) 
is  feasible  but,  at  some  subsequent  stage  z>l,  the  related  problem  (C‘)  beeomes 
infeasible.  The  reason  of  that  situation  is  that  we  eannot  determine  whether  the  souree  of 
infeasibility  is  that  the  (C)  is  aetually  integer-infeasible  (but  eontinuous-feasible),  or  that 
the  easeade  fixing  proeedure,  whieh  works  with  estimates  of  the  true  optimal  values  of 
the  variables,  makes  (C)  infeasible  unneeessarily.  We  ean  make  provisions  to  avoid  this 
situation  by  ereating  bigger  groups  of  tasks  (see  Eseudero  and  Salmeron  [2002]).  As  we 
will  show  in  Chapter  III,  in  the  test  eases  that  require  us  to  implement  F&R  ,  this  situation 
did  not  oeeur. 


A  similar  F&R  algorithm  for  our  (H)  formulation  has  been  implemented.  In  this 
ease,  the  partition  FJ  must  be  speeified  for  eaeh  goal  k  e  K  . 
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III.  IMPLEMENTATION  AND  COMPUTATIONAL  RESULTS 


In  this  chapter,  we  deseribe  some  details  of  our  implementation  and  present 
results  for  our  (H)  and  (C)  models  on  a  set  of  test  eases  drawn  from  Kirk  [1999]  and 
Arnold  [2000], 

A.  OVERVIEW 

We  implement  our  models  in  GAMS  [Brooke  et  al.  (1996)]  (version  2. 0.8. 3  with 
Revision  138  module)  using  the  CPLEX  solver  [GAMS-CPLEX  (2004)]  (version  9.0.0) 
on  a  Dell  Computer  Preeision  340  Pentium-4,  at  2.8  GHz  with  1  GB  of  random  aeeess 
memory. 

The  souree  eode  is  eurrently  divided  into  11  modules  eomprising  over  1,000  lines 
of  eode.  All  the  data  is  handled  separately  in  multiple  data  files.  Eor  speeifie  details  on 
the  use  of  the  eode,  see  TEAM  Eormulation  and  Code  [2004]  or  upeoming  updates  of  that 
doeument.  Next,  we  explain  some  of  the  details  speeified  through  our  data  fdes,  and 
show  examples  of  how  to  use  them. 

1.  Strategy  File 

The  “Strategy”  Data  File  (Figure  18)  allows  the  user  to  speeify  a  maximum 
solution  time,  the  number  of  F&R  stages  for  the  (C)  formulation,  and  whether  default 
weights  for  the  (C)  formulation  must  be  overwritten  by  eustomized  weights  (speeified 
separately)  or  not. 

The  default  weight  for  the  goal  in  the  k-th  position  on  the  goal  priority  list  is  1 0^'^. 
Sinee  there  are  eight  goals,  the  highest  priority  goal  on  the  list  will  reeeive  a  weight  of 
10^,  whereas  the  weight  for  the  lowest  priority  goal  is  1.  These  defaults  ean  be 
overwritten  by  the  user. 
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*  Parameters  for  the  strategy  for  CASE 


*  #  of  F&R  stages  when  a  weighed  objective  function  is  used 
Num_FR_stages  1 

*  Overwrite  default  weights 
Overwrite_weights  0 

*  Maximum  time  (seconds) 

Maxtime  600 


Figure  18.  Strategy  Data  File. 

This  file  eontains  user’s  strategy  information  to  solve  the  problem.  There  are  three  fields. 
The  first  two  (Num_FR_stages  and  Overwrite_weights)  apply  only  to  the  ease 
where  a  weighed  objeetive  funetion  is  used.  The  last  field  (Maxtime)  applies  to  both 
hierarehieal  and  eombined-goal  formulations. 


It  is  important  to  remark  that  the  use  of  adequate  goals  is  fundamental  to  aehieve 
the  desired  results.  Determining  weights  is  a  eritieal  task,  sinee  these  weights  must 
ensure  that  no  alloeation  ean  ever  produee  a  lower  priority  goal  to  outweigh  a  higher 
priority  goal.  Our  program  allows  the  user  to  manipulate  the  weights  in  order  to  avoid 
that  problem.  But,  the  validity  of  the  solutions  depends  on  the  extent  to  whieh  the 
operator  ean  determine  appropriate  weights  for  eaeh  problem.  This  involves  eonsidering 
possible  ranges  of  variation  for  all  goals  and  potential  sourees  of  numerieal  instability. 
These  are  two  examples  of  why  the  seleetion  of  weights  needs  to  be  earried  out  earefully: 

•  Suppose  the  first  goal  units  are  “number  of  unmet  task-parts”  and  the  seeond  goal 
units  are  “penalty  units  for  firing  unit  engagement  in  the  alloeation.”  Furthermore, 
assume  that  there  is  only  one  ship  available  and  one  task-part  to  be  alloeated.  If  the 
employment  penalty  of  the  ship  is  100,  the  default  weights  (wi=10  ,  W2=10  )  will  not 
eause  the  desired  effeet:  one  unmet  task-part  has  a  eost  of  10\l=10^,  whereas 
engaging  the  ship  in  the  alloeation  has  a  eost  of  10  x  100=  10  .  Consequently,  we 
would  eonelude  that  it  is  preferable  not  to  alloeate  the  task-part,  eontradieting  our 
initial  objeetive  of  giving  more  priority  to  the  alloeation  of  all  tasks.  To  avoid  this 
problem,  it  is  reeommended  to  normalize  the  objeetive  funetion  eoeffieients,  when 
possible,  for  example,  we  eould  have  used  an  employment  penalty  equal  to  one, 
instead  of  100.  In  other  eases,  it  may  be  neeessary  to  eustomize  all  of  the  weight 
values. 

•  Another  faetor  to  eonsider  is  the  optimality  toleranee  used  to  solve  the  models.  The 
default  value  is  0.01  (1%)  (this  ean  also  be  overwritten  in  by  the  user).  This  provokes 
an  early  termination  of  the  solving  proeess  (Braneh-and-Bound)  when  the  ineumbent 
solution  is  proven  to  be  within  the  speeified  toleranee  from  optimality.  Thus, 
objeetives  with  lower  weights  may  be  overlooked  beeause  ehanges  in  them  eannot 
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provoke  changes  in  the  objective  function  that  offset  the  tolerance  level.  An  example 
of  this  is  our  primary  validation  example  in  Section  II-D  (see  Figure  17).  There,  the 
lowest  priority  goal  optimal  value  is  10.0  (found  by  the  (H)  formulation).  However, 
when  this  goal  is  combined  with  the  others,  the  value  achieved  is  0.0.  In  practice,  this 
may  happen  even  if  the  optimality  tolerance  is  set  to  0%,  because  of  rounding  errors 
in  the  solver. 

The  Maxtime  field  is  the  allowed  time  to  solve  the  model  without  including  the 
generation  and  data  management  time,  which  may  not  be  negligible  in  some  cases. 
Unfortunately,  GAMS  suffers  from  enormous  overhead  for  model  generation,  because, 
for  example,  it  is  not  able  to  use  previously  generated  models  for  the  same  problem,  with 
perhaps  minor  changes  (such  as  changing  an  objective  function  coefficient,  or  adding  an 
additional  constraint).  Thus,  each  model  in  the  sequential  process  to  solve  the 
hierarchical  formulation  needs  to  be  regenerated,  and  the  same  occurs  with  each  model 
for  the  F&R  variable-partition  procedure.  While  we  realize  that,  at  this  moment, 
overhead  time  (as  part  of  the  total  time)  is  of  interest,  we  also  understand  that  much  of 
this  time  could  be  avoided  if,  for  example,  we  used  the  native  interfaces  for  the  different 
solvers  (e.g.,  CPLEX  callable  library,  OSL  EKK-callable  library,  etc.),  instead  of  using 
GAMS  to  create  the  models. 


Once  the  Maxtime  is  established,  the  following  time  allocation  is  made  for  the 
different  models: 

•  Combined  objective  function  with  one  E&R  stage: 

In  this  case,  all  of  Maxtime  seconds  are  allocated  to  solve  the  model. 


Hierarchical  objective  function  with  one  E&R  stage  per  goal: 


We  allocate  remaining  time  x 


0.5+- 


1 


#  of  remaining  goals  + 1 


seconds  to  the 


incumbent  goal.  Eor  example,  if  Maxtime  is  600  seconds  and  there  are  eight  goals,  then 
367  seconds  is  the  allotted  time  to  solve  the  first  goal.  If  it  can  be  solved  before  that  time 
limit,  for  example,  after  200  seconds,  then  there  are  400  seconds  left  for  the  remaining 
seven  goals.  Thus,  the  second  goal  receives  250  seconds,  and  so  forth. 
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Combined  objective  function  with  n  >  1  F&R  stages: 


( 


We  allocate  remaining  timex 


0.5+- 


1 


V 


seconds  to  the 


#  of  remaining  stages  +  \  j 
incumbent  stage.  The  idea  is  the  same  as  the  one  explained  in  the  above  bullet. 

•  Hierarchical  objective  function  with  nt  F&R  stages  for  goal  k,  where  >  1  for 


some  k  ^  K  . 

We  allocate  t^=  remaining  timex 
the  incumbent  goal,  k. 


f 


0.5+- 


1 


A 


#  of  remaining  goals  + 1 


seconds  to 


Then,  we  allocate  remaining  time^  x 


0.5+- 


1 


seconds 


#  of  remaining  stages + 1  y 
to  the  incumbent  stage  within  the  goal,  where  remaining  time,^  is  initialized  to  C  (total 


allotted  time  for  the  goal),  and  stages  are  for  goal  k. 

2.  Goal  Parameters  File 

Another  important  file  specifies  goal  data  (see  Figure  19): 


*  Goal  Parameters 


UNMET  TP 

Goal  Priority 

1 

Max  Deviation 
.20^ 

Num  FR  stages 

1 

Goal  weight 
1E7 

EMP  PEN 

2 

.20 

1 

1E5 

EXPEND  FU 

3 

.20 

1 

1E4 

LEVEL  FU 

4 

.20 

1 

1E3 

PR  SPREAD 

5 

.20 

1 

1E2 

BU  SPREAD 

6 

.20 

1 

lEl 

MCUBE  POS 

7 

.20 

1 

lEO 

SALVO  CAP 

8 

.20 

1 

lE-1 

Figure  19.  Goal  Parameters  Data  File. 

This  file  provides  flexibility  in  the  way  the  model  is  set  up.  There  are  four  columns  of 
data  for  each  goal.  The  first  field  (Goal_Priority)  applies  to  both  the  (C)  and  (H) 
formulations.  The  second  and  third  fields  (Max_Deviation  and  Num_FR_stages, 
respectively)  apply  to  the  hierarchical  formulation.  The  last  field  (Goal  weight) 
applies  to  the  combined-goal  formulation. 


The  Goal  Priority  field  must  contain  a  list  of  all  integers  from  one  through 
eight,  where  each  number  must  be  entered  exactly  once.  The  number  indicates  the 
priority  of  the  goal  and  will  be  used  for  the  following  purposes: 


46 


•  If  the  combined-goal  is  being  used,  a  default  weight  equal  to  10^'^  will  be  allocated 
for  the  goal  whose  priority  is  precisely  k.  These  weights  will  be  used  in  the  combined 
objective  function  unless  otherwise  specified  by  the  user  by  means  of  the 
Overwrite  weights  flag  (specified  in  the  aforementioned  Strategy  Data  File). 

•  If  a  hierarchical  formulation  is  being  used,  the  list  causes  the  highest  priority  goal 
(i.e.,  the  goal  whose  priority  is  equal  to  one)  to  be  optimized  first,  then  the  second 
(subject  to  a  restriction  on  the  first),  and  so  forth. 

The  Max_Deviation  field  for  a  goal  whose  priority  is  A:<8  specifies  the  level  of 
flexibility  in  the  achievement  required  for  this  goal  in  the  (H)  formulation.  It  sets  the 
maximum  deviation  from  the  goal’s  optimum  when  any  lower  priority  goal  k’  >  k  is 
being  optimized.  For  example,  using  the  data  in  Figure  19,  if  the  optimal  solution  for 
goal  #1  is  100,  then  goal  #2  will  be  optimized  subject  to  a  maximum  of  120  for  goal  #1. 
Remark:  The  deviation  must  be  specified  in  absolute  value  because  the  direction  for  each 
goal  (min.  or  max.)  is  already  implemented  in  the  code.  Thus,  if  goal  #1  were  a 
maximization  goal,  then  goal  #2  would  have  been  optimized  subject  to  a  minimum  of  80 
for  goal  #1. 

The  Num  FR  stages  field  contains  the  number  of  F&R  stages  for  each  goal.  If 
this  number  is  set  to  one,  the  optimization  problem  associated  with  that  goal  will  be 
solved  in  one  single  pass.  It  is  anticipated  that  some  goals  (such  as  “Spread  Primary 
Task-Parts”)  may  require  more  computational  work  than  others,  and  it  is  therefore 
advised  to  specify  more  than  one  F&R  stage  for  these  goals  in  large  cases. 

The  Goal  Weight  field  is  only  used  when  the  user  solves  the  combined-goal 
formulation  with  customized  weights.  The  Overwrite_weights  flag  (in  the 
Strategy.dat  file)  must  be  set  to  “1.”  Then,  the  Goal_Weight  field  can  be  used  to 
specify  new  weights  for  each  goal,  which  will  overwrite  the  default  ones. 

3,  Output  Files 

The  information  of  goal  achievement  and  time  used  depends  on  the  formulation 

used. 
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Total  solve  time: 

0 . 3  sec . 

G  1 

UNMET  TP 

(min) 

o 

o 

G  2 

EMP  PEN 

(min) 

1 . 0 

G  3 

EXPEND  FU 

(max) 

o 

o 

G  4 

LEVEL  FU 

(min) 

16.0 

G  5 

PR  SPREAD 

(max) 

1.0 

G  6 

BU  SPREAD 

(max) 

1 . 0 

G  7 

MCUBE  POS 

(min) 

25.0 

G  8 

SALVO  CAP 

(max) 

36.0 

OBJ_ 

VAL  W  : 

1159114 

o 

o 

OBJ_ 

EST  W  : 

1158614 

o 

o 

Figure  20.  Goal  Achievement  Results:  Combined-Goal  Formulation 


In  the  combined-goal  results,  the  output  (see  Figure  20)  shows: 

•  The  total  solve  time 

•  The  objective  achieved  for  each  goal 

•  The  combined-goal  achievement  and  a  lower  bound  on  the  best  possible  result. 
Remark:  the  combined-goal  formulation  is  stated  as  a  minimization  model,  see 
Section  II. B. 


Total  solve  time:  1.9  sec. 

MULTI-OBJECTIVE  TABLE 


G  1 

G  2 

G  3 

G  4 

G  5 

G  6 

G  7 

G  8 

Time 

(OBJ  VAL  H) (OBJ 

EST  H) 

(min) 

(min) 

(max) 

(min) 

(max) 

(max) 

(min) 

(max) 

(sec) 

G 

1: 

UNMET  TP 

0.0 

1.0 

0.0 

20.0 

0.0 

0.0 

29.0 

0.0 

0.2 

0.00 

0.00 

G 

2  : 

EMP  PEN 

0.0 

1.0 

0.0 

20.0 

0.0 

0.0 

29.0 

0.0 

0.2 

1.00 

1.00 

G 

3: 

EXPEND  FU 

0.0 

1.0 

0.0 

20.0 

0.0 

0.0 

29.0 

0.0 

0.3 

0.00 

0.00 

G 

4  : 

LEVEL  FU 

0.0 

1.0 

0.0 

16.0 

0.0 

0.0 

29.0 

0.0 

0.3 

16.00 

16.00 

G 

5: 

PR  SPREAD 

0.0 

1.0 

0.0 

18.0 

2.0 

0.0 

29.0 

0.0 

0.2 

2.00 

2.00 

G 

6: 

BU  SPREAD 

0.0 

1.0 

0.0 

18.0 

2.0 

0.0 

29.0 

0.0 

0.3 

0.00 

0.00 

G 

7  : 

MCUBE  POS 

0.0 

1.0 

0.0 

18.0 

2.0 

0.0 

25.0 

0.0 

0.2 

25.00 

25.00 

G 

8: 

SALVO  CAP 

0.0 

1.0 

0.0 

18.0 

2.0 

0.0 

25.0 

36.0 

0.2 

36.00 

36.00 

Figure  21 .  Goal  Achievement  Results:  Hierarchical  Formulation. 


When  the  (H)  formulation  is  used,  the  output  (see  Figure  21)  displays: 

•  The  total  solve  time 

•  The  objective  achieved  for  each  goal  at  each  level  of  the  hierarchical-optimization 
process 

•  The  time  needed  to  solve  each  goal 

•  A  bound  for  the  best  possible  value  for  each  goal 

For  example,  when  we  solve  for  the  first  goal  (minimization),  we  obtain  zero  for 
that  goal  and  other  non-optimal  values  for  the  other  goals  (see  the  line  headed  “G  1”). 
This  is  clearer,  for  example,  by  looking  at  goals  #3  and  number  #4.  When  goal  #3  is 
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optimized,  the  result  value  aehieved  for  goal  #3  is  zero  whereas  it  is  20  for  goal  #4. 
Then,  when  we  optimize  goal  #4,  we  are  able  to  reduee  this  value  to  16.  Moreover,  we 
know  this  value  is  optimal  beeause  the  bound  (last  eolumn  value)  indieates  that  the  best 
estimate  is  preeisely  16.  However,  when  we  optimize  goal  #5,  whieh  improves  its  value 
from  zero  to  two,  we  notiee  that  goal  #4  worsens  from  16  to  18.  This  is  due  to  the  faet 
that  we  have  allowed  a  maximum  of  20%  deviation  (see  Figure  19)  from  its  optimal  for 
this  goal. 

Of  eourse,  the  last  row  refers  to  the  last  goal  that  is  optimized  and  indieates  the 
values  aehieved  for  all  of  the  goals  after  the  optimization  proeess  is  eompleted.  As  we 
ean  see,  these  values  slightly  differ  from  those  aehieved  with  the  weighed  objeetive 
funetion,  due  to  the  degree  of  flexibility  allowed  in  some  goals.  If  an  optimal  solution  is 
required,  then  the  Max  Deviation  parameter  must  be  set  to  0.00  for  all  goals,  and  the 
optimality  toleranee  must  be  set  to  0.00  as  well. 

Figure  24  is  an  example  of  alloeation  output.  The  example  shows,  for  example, 
two  PR  task-parts  in  task  “Tl-1”  in  target  “Tl”  alloeated  to  two  III-C  weapons  in  half- 
modules  #1  and  #2  on  ship  DDG72.  Also,  we  ean  see  that  submarine  SSN775  is  engaged 
in  the  alloeation;  in  partieular,  two  II-C  weapons  and  one  III-C  weapon  from  the  torpedo 
room  are  used. 


G 

T 

TASKING  ALLOCATION 

PM  S 

H 

B 

T/C 

W 

Tl 

Tl-1 

PR 

TPMl 

DDG72 

h2 

III-C 

TPM2 

DDG72 

hi 

III-C 

RS 

TPMl 

DDG72 

h3 

III-C 

BU 

TPMl 

CG73 

h2 

III-C 

T2 

T2-1 

PR 

TPMl 

SSN755 

T 

II-C 

T2-2 

PR 

TPMl 

SSN755 

T 

III-C 

T2-3 

PR 

TPMl 

SSN755 

T 

II-C 

Figure  22.  Alloeation  Results 
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4,  Example 

The  flexibility  that  our  eode  allows  is  shown  in  the  following  example.  We 
eonsider  “Case  1”  drawn  from  Kirk’s  Thesis  [1999].  Initially,  we  eonsider  minimizing 
unmet  task-parts  as  the  highest  priority  goal,  and  maximizing  residual  salvo  as  the  least 
priority  goal.  However,  the  TSC  may  require  a  ehange  in  the  priorities,  and  he  or  she 
may  deem  maximizing  residual  salvo  as  the  seeond  priority  goal  (e.g.,  Kuykendall 
[1999]),  to  the  detriment  of  other  intermediate  goals.  We  eonduet  this  test  by  simply 
ehanging  the  priority  list  in  the  goal  data  fde,  as  explained  above.  As  we  show  in  Figure 
23,  the  aehievements  in  some  of  the  goals  ehange  in  these  two  different  eases. 


Goal  Name 

Direction 

Initial 

Goal 

Priority 

Value 

New  Goal 
Priority 

Value 

Unmet  Task-parts 

Min 

1 

0.0 

1 

0.0 

Employment  Penalty 

Min 

2 

0.0 

8 

1.0 

Expend  Firing  Units 

Max 

3 

N/A 

3 

N/A 

Level  Firing  Units 

Min 

4 

34.0 

4 

34.0 

Primary  Spread 

Max 

5 

3.0 

5 

4.0 

Back-up  Spread 

Max 

6 

N/A 

6 

N/A 

Min 

7 

92.0 

7 

92.0 

Residual  Salvo 

Max 

8 

159.0 

2 

160.0 

Figure  23.  Case  1:  Changing  Goal  Priorities. 


An  important  advantage  of  the  (H)  implementation  is  that  even  if  we  eannot  prove 
that  an  ineumbent  solution  is  optimal,  a  bound  on  the  best  possible  solution  is  provided. 
This  may  allow  the  TSC  to  let  more  solving  time  to  optimize  partieular  goal  in  a 
subsequent  run.  An  example  of  this  is  “Case  5  A”  from  Kirk  [1999].  As  we  ean  see  from 
Figure  24,  we  aehieve  38.5  for  Goal  #4  before  the  solving  time  limit  (set  to  400  seeonds 
for  this  test)  is  exeeeded.  At  this  point,  we  eannot  guarantee  that  38.5  is  the  optimal 
value  for  leveling  missiles  aeross  firing  units,  but  we  know  that  in  any  ease  we  eannot 
obtain  a  better  seore  than  37.4. 
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Goal 

Priority 

Goal  Name 

Direction 

(H) 

Formulation 

Best  Bound 

1 

Unmet  Task-parts 

Min 

0.0 

0.0 

2 

Employment  Penalty 

Min 

0.0 

0.0 

3 

Expend  Firing  Units 

Max 

72.0 

73.5 

4 

Level  Firing  Units 

Min 

38.5 

37.4 

5 

Primary  Spread 

Max 

4.0 

4.0 

6 

Back-up  Spread 

Max 

0.0 

0.0 

7 

M'’  List  Position 

Min 

594.0 

586.7 

8 

Residual  Salvo 

Max 

222.0 

242.8 

Figure  24.  Case  5A:  Comparison  of  Solution  Values  and  Best  Bounds. 

We  limit  the  solving  time  to  400  seconds  for  this  run.  Shaded  cells  shows  goals  for 
which  the  solution  achieved  has  not  been  proven  optimal. 


In  the  case  of  the  (C)  implementation,  we  do  not  have  a  measure  of  optimality  for 
each  individual  goal,  but  for  the  whole  combined  goal  that  makes  up  the  objective 
function. 

B,  SCENARIO  DESCRIPTION  AND  COMPUTATIONAL  RESULTS 

1,  Scenario  Description 

Our  test  focuses  on  existing  scenarios  developed  by  previous  NPS  students  and 
the  NSWCDD.  We  also  incorporate  two  new  scenarios  to  our  test  cases.  In  total,  we 
tested  23  cases,  including  sub-scenarios. 

The  first  six  scenarios  were  developed  by  NSWCDD  and  tested  initially  by  Kirk, 
and  then  by  Hodge  and  Arnold.  Some  of  these  test  scenarios  include  several  excursions, 
making  a  total  of  1 1  test  cases.  The  number  of  tasks  in  these  cases  ranges  from  13  tasks 
in  “Case  1”  to  104  tasks  in  “Case  6.”  All  of  Kirk’s  11  cases  only  include  ships  as  firing 
units,  and  do  not  have  multiple  missile  requirements  for  the  same  task-part,  among  other 
limitations.  For  additional  details,  consult  Kirk  [1999]. 

Arnold  adds  six  new  scenarios  to  test  new  features  such  as  submarines,  ghost 
tasks,  targets  containing  multiple  tasks,  etc.  For  details,  see  Arnold  [2000]. 

In  addition  to  these  twelve  scenarios  (and  their  variants),  we  create  two  more 
scenarios:  the  first,  “Case  13,”  is  strongly  based  on  Wingearf  s  most  complicated  case 
(Strike  Three  of  RIMPAC  00-3  Exercise),  which  has  been  recreated  manually  by  using 
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the  information  available  from  his  thesis  (as  opposed  to  Kirk’s  and  Arnold’s  scenarios  for 
which  electronic  data  were  available).  The  new  second  scenario,  “Case  14”  is  a 
hypothetical  case  intended  to  implement  multiple  specifications  (described  in  Chapter  II) 
in  a  single  case,  for  example,  submarine  back-up  tasking,  non-separable  tasking,  TSC 
specifications,  etc.  “Case  13”  includes  100  tasks  and  “Case  14”  has  26  tasks. 

2.  Test  Case  Set  Up 

We  test  all  scenarios  using  both  the  (C)  and  (H)  formulations.  Unless  otherwise 
specified,  we  use  the  following  set  of  parameters  to  run  all  test  cases: 

•  (H)  formulation  (multiple  submodels,  one  for  each  goal,  optimized  hierarchically): 

•  Optimality  tolerance:  0.01  (1%)  for  each  submodel. 

•  “Maxtime”  parameter:  (maximum  solve  time  for  the  overall 
process  including  all  submodels):  3,600  seconds 

•  “Num_FR”  parameter  (number  of  F&R  stages  per  submodel): 
One  (i.e.,  solve  the  submodel  as  a  whole  without  variable  partition) 
for  all  goals. 

•  “Max  Deviation”  parameter  (percentage  of  deviation  from  optimal 
goal  value):  0.01  (1%)  for  all  the  goals. 

•  (C)  formulation  (one  combined  objective  function  with  weighed  goals): 

•  Optimality  tolerance:  0.01  (1%). 

•  “Maxtime”  parameter:  3,600  seconds. 

•  “Num  FR”  parameter  (number  of  F&R  stages):  One  (i.e.,  solve 
the  model  as  a  whole  without  variable  partition). 

•  Weights  used:  =  10*^*,  Vk  e  &  =  {1,...,8} 

In  some  special  cases,  we  will  use  F&R  with  more  than  one  stage  in  order  to 
increase  the  model  or  submodel  solvability.  For  example,  we  may  be  applying  F&R  to 
the  submodel  associated  with  a  special  goal  in  the  (H)  formulation  or  to  the  model  in  the 
(C)  formulation. 

In  other  cases,  we  may  be  reducing  the  maximum  time  allotted  for  the  overall 
process  in  order  to  check  how  this  affects  the  quality  of  the  incumbent  solution  as  well  as 
the  best  bound  on  the  optimal  solution  if  this  is  not  achieved. 
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3,  Results 

We  compare  our  results  to  those  obtained  by  Kirk  [1999]  and  Arnold  [2000]  (the 
former  for  applicable  cases  only).  We  omit  Hodge  [1999]  from  the  comparasion  because 
his  results  are  always  outperformed  by  Kirk  and  Arnold. 

Table  1  shows  the  results  for  all  scenarios.  The  shaded  cells  show  the  key 
differences  between  our  solutions  and  the  existing  ones.  As  explained  earlier  in  this 
thesis,  we  display  solution  time  as  that  required  to  solve  the  model,  omitting  overhead 
time  caused  by  model  generation  and  other  data  management  in  GAMS. 

For  the  first  six  scenarios  and  their  variants,  we  notice  the  reduction  in 
computational  time  to  solve  our  (H)  and  (C)  formulations  when  compared  to  that  reported 
by  Kirk  [1999],  who  produced  the  last  available  MIP  formulation  prior  to  this  work.  One 
of  the  most  dramatic  cases  is  “Case  5,”  in  which  solution  time  decreases  from  10,976 
seconds  for  Kirk’s  Hierarchical  solution  to  40.5  seconds  for  our  (H)  formulation. 

Regarding  the  quality  of  the  solutions  provided  by  our  formulations,  it  is  the  same 
as  or  better  than  that  of  Kirk  and  Hodge  (where  available),  and  Arnold  (in  all  cases, 
except  “Case  2”  where  we  believe  there  is  an  error  in  the  result  produced  by  Arnold, 
which  is  also  confirmed  by  Kirk’s  result).  (Remark:  In  some  cases,  such  as  “Case  IB” 
using  the  (H)  formulation,  we  notice  that  the  optimal  value  reported  for  Goal  #8  is  159, 
whereas  the  optimal  is  160.  This  and  other  minor  discrepancies  can  be  attributed  to  the 
optimality  tolerance  established  at  1%  for  all  models.  In  these  cases,  we  still  consider 
both  solutions  of  the  same  quality.) 

“Case  13”  is  similar  to  that  of  LT  Wingeart’s  most  complicated  case:  RIMPAC 
00-3  Strike  Three  case,  with  100  tasks  [Wingeart,  2001].  Because  the  full  data  set  for 
this  case  was  not  available  at  the  time  of  developing  this  work,  we  manually  recreated  the 
case  using  the  information  shown  in  Wingerat’s  thesis.  Consequently,  our  results  can  be 
considered  comparable  to  those  of  Wingeart  in  terms  of  scope  and  case  difficulty,  but  not 
necessarily  on  a  goal  by  goal  basis. 

Finally,  we  tested  a  completely  new  scenario,  “Case  14,”  to  further  validate  the 
incorporation  of  all  problem  specifications  into  our  formulation. 
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Test 

Case 

Implementation 

Source 

GOALS 

Solution 

Time 

(Seconds) 

#  1 

(Min): 

Unassigned 

Task-parts 

#2 

(Min): 

Penalty 

Units* 

#3 

(Max): 

Expend 

Units 

#4 

(Min): 

Deviation 

Level 

#5 

(Max): 

Primary 

Missiles 

#6 

(Max): 

Back-up 

Missiles 

#7 

(Min): 

M^  List 
Position 

#8 

(Max): 

Residual 

Salvo 

1 

Kirk  (Hier.  R.) 

0 

0 

N/A 

34 

3 

N/A 

92 

160 

26 

Arnold  (Heur.) 

0 

0 

N/A 

34 

3 

N/A 

92 

160 

2 

(C)  Model 

0 

0 

N/A 

34 

3 

N/A 

92 

160 

0.3 

(H)  Model 

0 

0 

N/A 

34 

3 

N/A 

92 

159 

1.8 

lA 

Kirk  (Hier.  R.) 

0 

N/A 

10 

8.7 

N/A 

N/A 

92 

160 

34 

Arnold  (Heur.) 

0 

N/A 

10 

8.7 

N/A 

N/A 

92 

160 

2 

(C)  Model 

0 

N/A 

10 

8.7 

N/A 

N/A 

92 

160 

4.3 

(H)  Model 

0 

N/A 

10 

8.7 

N/A 

N/A 

92 

160 

1.8 

IB 

Kirk  (Hier.  R.) 

0 

N/A 

10 

8.7 

N/A 

2 

92 

160 

36 

Arnold  (Heur.) 

0 

N/A 

10 

8.7 

N/A 

3 

92 

160 

2 

(C)  Model 

0 

N/A 

10 

8.7 

N/A 

4 

92 

160 

0.3 

(H)  Model 

0 

N/A 

10 

8.7 

N/A 

4 

92 

159 

2.1 

2 

Kirk  (Hier.  R.) 

0 

0 

N/A 

23 

3 

N/A 

184 

157 

335 

Arnold  (Heur.) 

0 

0 

N/A 

13.5** 

2 

N/A 

221 

162 

3 

(C)  Model 

0 

0 

N/A 

23 

3 

N/A 

184 

157 

0.8 

(H)  Model 

0 

0 

N/A 

23 

3 

N/A 

184 

157 

5.6 

2A 

Kirk  (Hier.  R.) 

0 

1 

N/A 

42.5 

2 

N/A 

196 

146 

103 

Arnold  (Heur.) 

0 

1 

N/A 

44.5 

4 

N/A 

196 

141 

4 

(C)  Model 

0 

1 

N/A 

42.5 

4 

N/A 

186 

146 

0.4 

(H)  Model 

0 

1 

N/A 

42.5 

4 

N/A 

186 

145 

1.9 

3 

Kirk  (Hier.  R.) 

0 

0 

N/A 

15.5 

3 

N/A 

368 

129 

1625 

Arnold  (Heur.) 

0 

1 

N/A 

1.5 

3 

N/A 

368 

130 

5 

(C)  Model 

0 

0 

N/A 

15.5 

3 

N/A 

368 

140 

1.5 

(H)  Model 

0 

0 

N/A 

15.5 

3 

N/A 

368 

140 

14.5 

4 

Kirk  (Hier.  R.) 

0 

1 

N/A 

0 

4 

N/A 

496 

127 

6133 

Arnold  (Heur.) 

0 

1 

N/A 

16 

4 

N/A 

508 

107 

10 

(C)  Model 

0 

1 

N/A 

0 

4 

N/A 

484 

127 

31.7 

(H)  Model 

0 

1 

N/A 

0 

4 

N/A 

488 

126 

23.1 

5 

Kirk  (Hier.  R.) 

0 

0 

N/A 

11.4 

4 

N/A 

558 

268 

10976 

Arnold  (Heur.) 

0 

0 

N/A 

25.4 

6 

N/A 

599 

252 

10 

(C)  Model 

0 

0 

N/A 

11.4 

4 

N/A 

558 

268 

6.1 

(H)  Model 

0 

0 

N/A 

11.4 

4 

N/A 

562 

268 

40.4 

5A 

Kirk  (Hier.  R.) 

0 

0 

65 

34 

5 

N/A 

594 

242 

2356 

Arnold  (Heur.) 

0 

0 

73 

40 

4 

N/A 

596 

227 

10 

(C)  Model 

0 

0 

73 

40 

4 

N/A 

597 

230 

4.5 

(H)  Model 

0 

0 

72 

38.5 

4 

N/A 

597 

241 

3367.9 

5B 

Kirk  (Hier.  R.) 

0 

1 

N/A 

21.4 

4 

N/A 

564 

268 

34340 

Arnold  (Heur.) 

0 

1 

N/A 

35.1 

5 

N/A 

629 

240 

16 

(C)  Model 

0 

1 

N/A 

21.4 

5 

N/A 

583 

268 

1349 
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Test 

Case 

Implementation 

Source 

GOALS 

Solution 

Time 

(Seconds) 

#  1 

(Min): 

Unassigned 

Task-parts 

#2 

(Min): 

Penalty 

Units* 

#3 

(Max): 

Expend 

Units 

#4 

(Min): 

Deviation 

Level 

#5 

(Max): 

Primary 

Missiles 

#6 

(Max): 

Back-up 

Missiles 

#7 

(Min): 

M^  List 
Position 

#8 

(Max): 

Residual 

Salvo 

(H)  Model 

0 

1 

N/A 

20.6 

5 

N/A 

592 

270 

2328 

6 

Kirk  (Hier.  R.) 

0 

0 

69 

24 

6 

N/A 

904 

209 

261 

Arnold  (Heur.) 

0 

0 

69 

24 

6 

N/A 

902 

201 

7 

(C)  Model 

0 

0 

78 

13.5 

5 

N/A 

892 

240 

22.2 

(H)  Model 

0 

0 

78 

13.5 

5 

N/A 

896 

238 

226.8 

7 

Arnold  (Heur.) 

0 

0 

N/A 

34 

2 

N/A 

92 

160 

2 

(C)  Model 

0 

0 

N/A 

34 

3 

N/A 

92 

152 

0.3 

(H)  Model 

0 

0 

N/A 

34 

3 

N/A 

92 

152 

1.4 

7A 

Arnold  (Heur.) 

0 

0 

10 

34.7 

N/A 

N/A 

104 

152 

5 

(C)  Model 

0 

0 

10 

34.7 

N/A 

N/A 

107 

147 

1.7 

(H)  Model 

0 

0 

10 

34.7 

N/A 

N/A 

107 

147 

1.4 

7B 

Arnold  (Heur.) 

0 

0 

10 

8.7 

N/A 

3 

92 

160 

2 

(C)  Model 

0 

0 

10 

8.7 

N/A 

4 

92 

150 

6 

(H)  Model 

0 

0 

10 

8.7 

N/A 

4 

92 

150 

1.9 

8 

Arnold  (Heur.) 

0 

0 

N/A 

36 

N/A 

N/A 

80 

162 

2 

(C)  Model 

0 

0 

N/A 

36 

N/A 

N/A 

80 

162 

0.3 

(H)  Model 

0 

0 

N/A 

36 

N/A 

N/A 

80 

161 

2 

8A 

Arnold  (Heur.) 

0 

0 

N/A 

36 

N/A 

N/A 

81 

162 

2 

(C)  Model 

0 

0 

N/A 

36 

N/A 

N/A 

80 

153 

0.2 

(H)  Model 

0 

0 

N/A 

36 

N/A 

N/A 

80 

153 

1.4 

9 

Arnold  (Heur.) 

0 

0 

N/A 

34 

3 

N/A 

92 

160 

3 

(C)  Model 

0 

0 

N/A 

34 

3 

N/A 

92 

160 

0.3 

(H)  Model 

0 

0 

N/A 

34 

3 

N/A 

92 

160 

1.8 

9A 

Arnold  (Heur.) 

0 

0 

N/A 

34 

N/A 

3 

92 

160 

3 

(C)  Model 

0 

0 

N/A 

34 

N/A 

3 

92 

153 

0.3 

(H)  Model 

0 

0 

N/A 

34 

N/A 

3 

92 

153 

1.6 

10 

Arnold  (Heur.) 

0 

0 

N/A 

39 

3 

N/A 

157 

147 

3 

(C)  Model 

0 

0 

N/A 

31 

3 

N/A 

149 

152 

0.6 

(H)  Model 

0 

0 

N/A 

31 

3 

N/A 

149 

152 

2.9 

lOA 

Arnold  (Heur.) 

0 

0 

N/A 

37 

3 

N/A 

167 

147 

3 

(C)  Model 

0 

0 

N/A 

31 

3 

N/A 

149 

152 

0.6 

(H)  Model 

0 

0 

N/A 

31 

3 

N/A 

149 

152 

2.6 

11 

Arnold  (Heur.) 

0 

N/A 

0 

1 

N/A 

58 

162 

1 

(C)  Model 

0 

0 

N/A 

34 

3 

N/A 

92 

160 

0.3 

(H)  Model 

0 

0 

N/A 

34 

3 

N/A 

92 

159 

1.7 

12 

Arnold  (Heur.) 

0 

0 

N/A 

43 

7 

N/A 

654 

256 

12 

(C)  Model 

0 

0 

N/A 

15.8 

6 

N/A 

635 

272 

12 

(H)  Model 

0 

0 

N/A 

15.8 

6 

N/A 

641 

272 

46.4 

13 

Wing.(Heur(B)) 

0 

N/A 

N/A 

49 

8 

6 

128 

168 

N/A 

Wing.(Heur(W)) 

4 

N/A 

N/A 

73.5 

8 

3 

111 

177 

N/A 
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Test 

Case 

Implementation 

Source 

GOALS 

Solution 

Time 

(Seconds) 

#  1 

(Min): 

Unassigned 

Task-parts 

#2 

(Min): 

Penalty 

Units* 

#3 

(Max): 

Expend 

Units 

#4 

(Min): 

Deviation 

Level 

#5 

(Max): 

Primary 

Missiles 

#6 

(Max): 

Back-up 

Missiles 

#7 

(Min): 

M^  List 
Position 

#8 

(Max): 

Residual 

Salvo 

Wing.(Split(B)) 

0 

N/A 

N/A 

49 

8 

5 

128 

169 

N/A 

Wing.(Split(W)) 

6 

N/A 

N/A 

73.5 

8 

3 

112 

176 

N/A 

(C)  Model 

0 

N/A 

N/A 

45 

6 

7 

448 

371 

7.3 

(H)  Model 

0 

N/A 

N/A 

45 

6 

7 

451 

371 

199 

14 

(C)  Model 

0 

0 

24 

7.3 

2 

2 

227 

178 

0.5 

(H)  Model 

0 

0 

24 

7.3 

2 

2 

227 

177 

2.9 

Table  1.  Result  List. 

Maximum  time  allowed  is  3,600  seeonds  and  optimality  toleranee  is  0.01  for  all  runs. 

*;  To  be  eonsistent  with  reports  by  previous  students,  aetual  employment  penalties  are 
divided  by  100. 

**:  Arnold  solution  for  this  ease  appears  to  be  ineorreet. 

***:  Available  data  for  this  ease  appears  to  be  ineonsistent  with  Arnold’s  report. 

4.  Refinements  to  Tackle  the  Hardest  Test  Cases 

While  MIP  solution  quality  is  expeeted  to  exeeed  heuristies,  there  are  a  few 
instanees  in  whieh  solving  the  MIP  may  require  a  prohibitive  eomputational  time.  It  is 
not  possible  to  know  whieh  eases  will  require  an  unaeeeptable  time  to  be  optimally 
solved,  but  we  have  made  provisions  to  speed  up  the  MIP  solver. 

First,  we  establish  priority  branehing  as  follows:  Half-module-related  variables 
are  branehed  first.  Then,  torpedo-tube-related  variables  are  branehed,  and  lastly, 
eanister-eapsule -variables.  Within  eaeh  of  these  eategories,  the  more  weapons  of  a  given 
type  the  asset  has,  the  higher  the  branehing  priority.  For  example,  given  as  the 
number  of  weapons  of  type  w  in  half-module  h  on  ship  s,  the  priority  value  given  to 
^Iwtpm  is  +  Remark:  The  lower  the  priority  value,  the  sooner  the  assoeiated 

variable  is  branehed.  By  using  these  priorities,  we  notiee  that  a  slight  improvement  in 
eomputational  time  for  some  eases,  (typieally  under  10%). 
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In  this  section  we  report  the  results  obtained  by  (a)  reducing  the  total  solving  time 
and  (b)  using  F&R  with  several  stages  in  order  to  partition  our  models  into  easier 
subproblems. 

a.  Reducing  the  Solving  Time 

We  observe  that  “Case  5B”  is  one  of  the  most  difficult  cases  to  solve.  In 
fact,  when  we  let  the  maximum  solving  time  increase,  this  case  uses  up  all  the  time 
available.  When  we  decrease  the  maximum  solving  time  to  600  seconds,  the  (H) 
formulation  gives  the  same  near-optimal  solution  as  before  (see  Figure  25).  We  notice 
that  in  both  cases,  the  best  available  bound  for  goal  #4  is  11 .4,  while  the  best  solution 
found  is  20.6.  It  is  clear  that  this  particular  problem  is  consuming  all  the  available  time 
allocated  for  that  goal,  which  is  dynamically  computed  as  explained  in  Section  A. 


Goal 

Maximum  Time 
(36000  seconds) 

Priority 

Name 

Direction 

(H) 

Formulation 

Best 

Bound 

(H) 

Formulation 

Best 

Bound 

1 

Unmet  Task-parts 

Min 

0.0 

0.0 

0.0 

0.0 

2 

Employment  Penalty 

Min 

1.0 

1.0 

1.0 

1.0 

3 

Expend  Firing  Units 

Max 

N/A 

N/A 

N/A 

N/A 

4 

Level  Firing  Units 

Min 

20.6 

11.4 

20.6 

11.4 

5 

Primary  Spread 

Max 

5.0 

5.0 

5.0 

5.0 

6 

Back-up  Spread 

Max 

N/A 

N/A 

N/A 

N/A 

7 

Min 

592.0 

582.4 

592.0 

582.4 

8 

Residual  Salvo 

Max 

270.0 

271.0 

270.0 

271.0 

Solution  Time  (Seconds) 

21767.0 

584.4 

Figure  25.  Case  5B:  Reducing  Maximum  Time. 

We  set  the  maximum  solving  time  to  36,000  seconds  and  solve  the  case  with  the  (H) 
formulation.  The  solution  and  bounds  achieved  are  the  same  as  those  achieved  when  the 
maximum  time  is  set  to  600  seconds.  Some  goals  such  as  goal  #4  are  using  all  the  time 
allocated  to  solve  them. 

b.  Using  Fix  and  Relax 

In  Figure  26,  we  show  the  results  obtained  by  applying  F&R  to  our  (H) 
formulation  for  “Case  5  A.”  Without  F&R  ,  the  (H)  formulation  consumes  almost  all  the 
available  time  (3,600  seconds).  As  in  “Case  5B”  above,  some  of  the  goals  are  using  up 
all  the  allotted  time.  When  we  use  F&R  for  selected  goals  (identified  as  complicated 
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goals  for  this  case),  we  observe  that  the  total  solving  time  is  redueed  drastieally,  without 
eompromising  the  quality  of  the  solution  value  and  the  assoeiated  bounds  in  this  ease. 


Goal 

(H)  Formulation 

(H)  Formulation  with  F&R 

Priority 

Name 

Direction 

Value 

Best 

Bound 

H 

Value 

Best 

Bound 

1 

Unmet  Task-parts 

Min 

0.0 

0.0 

1 

0.0 

0.0 

2 

Employment  Penalty 

Min 

0.0 

0.0 

1 

0.0 

0.0 

3 

Expend  Firing  Units 

Max 

72.0 

73.5 

4 

73.0 

73.0 

4 

Level  Firing  Units 

Min 

38.5 

37.4 

4 

40.0 

40.0 

5 

Primary  Spread 

Max 

4.0 

4.0 

4 

4.0 

4.0 

6 

Back-up  Spread 

Max 

N/A 

N/A 

4 

N/A 

N/A 

7 

Min 

597.0 

586.7 

4 

597.0 

592.6 

8 

Residual  Salvo 

Max 

241.0 

242.4 

4 

239.0 

239.8 

Solution  Time  (Seconds) 

3367.9 

316.7 

Case  5A;  Results  with  Fix-and-Relax. 


Solution  quality  is  not  eompromised  in  this  ease. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER 

RESEARCH 


This  thesis  consolidates  the  latest  specifications  for  TLAM  pre-designation  into  a 
revised  mathematical  formulation,  and  implements  it  as  a  multi-objective  mixed- integer 
optimization  problem.  We  show  that  our  models  can  be  optimally  solved  in  reasonable 
time  for  most  of  the  test  cases,  and  in  those  cases  where  computational  time  is  an  issue, 
we  establish  mechanisms  to  alleviate  the  model  complexity  and  achieve  solutions  of 
acceptable  quality. 

On  the  other  hand,  through  our  models  we  validate  heuristic  results  for  a  number 
of  test  cases  developed  by  former  NFS  students.  Although  the  heuristic  may  fail  to 
provide  an  optimal  answer  even  for  the  highest  priority  goal,  its  overall  performance  can 
be  considered  satisfactory. 

There  are  several  aspects  in  which  this  thesis  contributes  to  and  enhances  earlier 
work  on  automated  TLAM  predesignation.  With  respect  to  Kirk’s  work  [Kirk,  1999],  the 
contributions  lie  on; 

•  Reformulating  multiple  parts  of  the  model,  making  it  possible  to  attain  the  solution  of 
problems  that  could  not  be  solved  before  by  the  combined-goal  formulation  (C),  and 
drastically  reducing  the  time  for  the  hierarchical  formulation  (H). 

•  Incorporating  all  new  specifications,  such  as  submarines,  targets,  multiple  task-part 
requirements,  etc. 

•  Implementing  a  Fix-and-Relax  (F&R)  variable-partition  procedure  to  help  solve  the 
most  complicated  cases. 

•  Implementing  our  programs  in  a  much  more  amenable  fashion,  where  the  front-end 
for  the  user  is  composed  of  a  number  of  data  files  completely  independent  of  the 
source  code.  For  example,  changing  goal  priorities  is  as  simple  as  updating  an  array 
of  eight  numbers.  All  of  our  test  cases  use  the  same  code,  and  all  reasonable  data  and 
option  customization  is  done  at  the  level  of  data  files. 

With  respect  to  Arnold’s  work,  our  (C)  and  (H)  formulations: 

•  Enhance  the  quality  of  the  solution  obtained. 

•  Provide  a  valid  bound  on  the  optimal  solution  in  case  this  is  not  achieved. 

•  Have  the  flexibility  to  arrange  goals  in  any  desired  order,  to  surrogate  certain  goals  to 
desired  levels  of  achievement  in  other  goals  ((H)  formulation),  or  to  allocate  different 
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weights  in  order  to  generate  a  variety  of  effieient  points  for  the  multi-objeetive 
problem  ((C)  formulation).  These  are  important  limitations  in  the  heuristie. 

The  work  in  this  thesis  is  part  of  an  ongoing  effort  by  students  and  faeulty  at  NFS 
and  NSWCDD.  In  that  eontext,  we  believe  that  further  validation  of  our  models  with  real 
data  (not  available  at  the  time  of  eompleting  this  thesis)  is  neeessary. 

We  also  realize  that  GAMS  may  not  be  the  most  appropriate  environment  for 
model  generation.  In  some  eases,  the  model  generation  and  exeeution  time  is  even 
greater  than  the  time  to  solve  the  (already  generated)  model  in  CPLEX.  This  burden 
intensifies  when  multiple  models  (with  minor  differenees  between  them)  need  to  be 
solved,  sueh  as  in  the  (H)  formulation  and/or  when  F&R  is  used. 
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APPENDIX  A.  FULL-TLAM  ALLOCATION  MODEL 


NOTATION 
1.  Indices 

H, 

B, 

G, 

T, 

P, 

M, 

W, 

I, 

A, 

K, 


2. 


Sets 


rp(j 

t!. 


Gf 

G 

^PR 


a  ’ 
R 


i^Exp 

nNS 


W 


set  of  ships;  s 

set  of  half-modules;  /?  e  //  =  |l,2,3,...} 

set  of  submarines;  h 
set  of  targets;  g  e  G 
set  of  tasks;  t  e  T 

set  of  task-parts;  p  =  [PR,RS,  BU,  GHj 

set  of  task-part  missile  requirements;  m,m'  e  M  =  |l,2,3,...} 

set  of  weapon  types;  weW 

set  of  subsets  near-simultaneous  tasks;  z  e  / 

set  of  launch  areas;  a  e  A 

set  of  goals;  k  ^  K  =  |l,2,...,8} 


subset  of  tasks  for  target  g;  cz  T 

subset  of  near-simultaneous  tasks  in  group  z;  T‘  ex  T 

subset  of  targets  which  are  geo-feasible  from  area  a;  G^  a  G 

subset  of  targets  which  are  designated  as  restricted;  G^  (^G 

subset  of  areas  whose  primary  task-parts  are  desired  to  be 
assigned  to  as  many  firing  units  as  possible;  ^4™  c  A 

subset  of  areas  whose  back-up  task-parts  are  desired  to  be 
assigned  to  as  many  firing  units  as  possible;  A^’^  cz  A 

subset  of  expend  ships;  5*^^^  c  S 
subset  of  expend  submarines;  a  B 

subset  of  submarines  that  are  desired  to  be  assigned  to  back-up 
task-part  allocation  of  any  task  t;  B^^  cz  B 

subset  of  task-parts  in  task  t  that  are  not  separable;  P,^^  cz  P 

subset  of  weapons  that  can  be  used  for  task  t;  Wj  cz  W 
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3. 


TSC^^, 


Data 


NExp 

fl  , 


subset  of  six-tuples  (s,h,w,t,p,m)  sueh  that  the  TSC  speeifies 
that  task-  part  missile  {t,p,m)  must  be  alloeated  to  weapon  type 
w  in  half-module  h  on  ship  s. 

subset  of  five-tuples  {b,w,t,p,m)  sueh  that  the  TSC  speeifies 
that  task-  part  missile  {t,p,m)  must  be  alloeated  to  weapon  type 
w  in  the  torpedo  room  of  submarine  b. 

subset  of  five-tuples  {b,w,t,p,m)  sueh  that  the  TSC  speeifies 
that  task-part  missile  {t,p,m)  must  be  alloeated  to  weapon  type 
w  in  the  eanister  eapsules  of  submarine  b. 

number  of  weapons  of  type  w  in  half-module  b  on  ship  5 

number  of  weapons  of  type  w  in  the  torpedo  room  of 
submarine  b 

number  of  weapons  of  type  w  in  the  eanister  eapsules  of 
submarine  b 

number  of  half-modules  on  ship  s 
number  of  torpedo  tubes  on  submarine  b 
number  of  missiles  required  for  task-part  p  of  task  t 
employment  penalty  for  ship  s 
employment  penalty  for  submarine  b 

number  of  non-expend  firing  units:  =  js*  \  + 15  \  | 

area  for  ship  5 
area  for  submarine  b 

weight  of  task-part  p  to  mateh  the  M  list  position  of  for  task- 
part  p 

3  T 

M  position  of  weapon  for  task  t 

relative  value  of  weapon  w  for  residual  salvo  eomputation 

maximum  number  of  firing  units  to  whieh  a  restrieted  target  g 
ean  be  alloeated 

pereentage  of  deviation  (from  the  best  possible)  allowed  for 
goal  k  when  a  less  priority  goal  is  eonsidered  in  the  objeetive 
funetion 
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direction  of  goal  k\  =  1  if  maximization,  /(-=-!  if 
minimization 

weight  given  to  goal  k  when  a  combined  goal  objective  is 
considered,  where  Oj  »  (»2  ^  ^  • 

4.  Derived  Data  and  Sets 


5. 


s:. 


B. 


target  for  task  t  (i.e.,  t  eT^) 
subset  of  geo-feasible  ships  for  task  t;  Sf  =  1 5 


‘  ’  subset  of  geo-feasible  submarines  for 
t;  =  |b  |t  e  For  a  =  ' 

Decision  Variables 


task 


^  shwtpm  ’ 


X^T 
bwtpm  ’ 


X^c 
^  bwtpm  ’ 


jS 
stpm  ’ 


Y 

btpm  ’ 


U. 


tpm  ’ 


ZZ" , 

sg’ 


zz 


bg  ’ 


E 


b  ’ 
5 


equals  1  if  a  weapon  of  type  w  in  half-module  h  of  ship  5  is 
used  to  allocate  missile  m  of  task-part  p  in  task  t,  0  otherwise 

equals  1  if  a  weapon  of  type  w  in  the  torpedo  room  of 
submarine  b  is  used  to  allocate  missile  m  of  task-part  p  in  task 
t,  0  otherwise 

equals  1  if  a  weapon  of  type  w  in  the  canister  capsules  of 
submarine  b  is  used  to  allocate  missile  m  of  task-part  p  in  task 
t,  0  otherwise 

equals  1  if  missile  m  of  task-part  p  in  task  t  is  allocated  to  ship 
5,  0  otherwise 

equals  1  if  missile  m  of  task-part  p  in  task  t  is  allocated  to 
submarine  b,  0  otherwise 

equals  1  if  task  t  is  unmet,  0  otherwise 

equals  1  if  ship  5  is  engaged  in  any  tasking  in  target  g,  0 
otherwise 

equals  1  if  ship  5  is  used,  0  otherwise 

equals  1  if  submarine  b  is  engaged  in  any  tasking  in  target  g,  0 
otherwise 

equals  1  if  submarine  b  is  used,  0  otherwise 

total  number  of  weapons  allocated  to  primary  task-parts  from 
expend  ships 

total  number  of  weapons  allocated  to  primary  task-parts  from 
expend  submarines 
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NS^ 

stp 

9 

RNl 

RNl 

RAVG, 

RDs, 

RDl 

prU 

prL 

BU^ 

BU!, 


D. 


shw  ’ 


D, 


BT 
hw  ’ 


D. 


BC 
bw  ’ 


equals  1  if  all  missile  requirements  for  non-separable  task-part 
p  of  task  t  are  assigned  to  ship  5,  0  otherwise 

equals  1  if  all  missile  requirements  for  non-separable  task-part 
p  of  task  t  are  assigned  to  submarine  b,  0  otherwise 

number  of  missiles  on  non-expend  ship  s  that  are  not  alloeated 
to  a  primary  task-part 

number  of  missiles  on  non-expend  submarine  b  that  are 
alloeated  to  a  primary  task-part 

average  of  the  residual  number  of  missiles  aeross  non-expend 
firing  units  after  alloeation 

differenee  between  the  residual  number  of  missiles  on  ship  s 
2inARAVG 

differenee  between  the  residual  number  of  missiles  on 
submarine  b  and  RA  VG 

equals  1  if  ship  5  is  assigned  to  alloeated  any  primary  task-part, 
0  otherwise 

equals  1  if  submarine  b  is  assigned  to  alloeated  any  primary 
task-part,  0  otherwise 

equals  1  if  ship  5  is  assigned  to  any  baek-up  task-part,  0 
otherwise 

equals  1  if  submarine  b  is  assigned  to  any  baek-up  task-part,  0 
otherwise 

equals  1  if  one  or  more  weapons  of  type  w  remain  on  ship  s  in 
half-module  b  after  firing  all  primary  task-part  missiles,  0 
otherwise 

number  of  weapons  of  type  w  remain  in  the  torpedo  room  of 
submarine  b  after  firing  all  primary  task-part  missiles,  0 
otherwise 

number  of  weapons  of  type  w  remain  in  the  eanister  eapsules 
of  submarine  b  after  firing  all  primary  task-part  missiles,  0 
otherwise 

objeetive  funetion  value  for  goal  k 
eombined-goal  objeetive 
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B, 


OBJECTIVE  FUNCTION 

Our  multi-objective  function  consists  of  the  following  goals: 


Goal  k=  \  :  Minimize  the  total  number  of  unmet  task-parts: 

Z;  =  minZi,  where  Zj  ^  J 

t  p  m=\ 

Goal  k  =  2:  Minimize  the  use  of  penalty  firing  units: 

Z;  =  mmZ2,  where  Z2  = 

s  b 

Goal  k  =  3:  Maximize  the  number  of  weapons  allocated  from  designated  expend  firing 
units: 

Zj  =  maxZ3,  where  Z3  = 

Goal  k  =  4:  Minimize  the  deviation  from  average  of  residual  number  of  missiles  on  non- 
expend  firing  units: 

Zj  =  minZ4,  where  Z4  =  ^  ^  RD^ 


Goal  k  =  5:  Maximize  the  number  of  firing  units  that  have  been  assigned  to  allocate 
primary  task-parts  in  launch  areas  so  designated: 


Z*  =  maxZj,  where  Z^ 


E 

a^PR'^ 


=a 


b\au  =a 


Goal  k=  6:  Maximize  the  number  of  firing  units  that  have  been  assigned  to  allocate  back¬ 
up  task-parts  in  launch  areas  so  designated: 


Z*  =  maxZg,  where  Zg  =  Y 


aeBU-* 


E  +  E 


h\aj^  =a 


Goal  k  =  1:  Minimize  the  missile  to  mission  matching  (M^)  list  position  of  all  selected 
weapons: 
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mf 

Z;  =minZ,,  where  ^7=ZZZ  ZZZ 

s  h=\  w|n^>0  JseS/’,  P  m=l 

\w^wj 

M 

Z  Z  Z  Z  Z  K^i^Zpm 

b  w|nfj>0  \b^Bj ,  P  m=l 

W£»'/ 

M 

Z  Z  Z  Z  Z 

b  w\nZ>0  \beBj' ,  P  >n=l 


Goal  A:  =  8:  Maximize  residual  salvo  capability: 


Z;  =  max Zg,  where  Z,  =  ^  ^  v^vAL  +  Z  Z  ^wAT+Z  Z 

S  A=1  M'|nfj„>0  b  w|nfj>0  *  w|nfir>0 


Combined  goal:  Minimize  a  weighed  average  of  the  goals  k  ^  K : 


Z;  =  minZo,  where  Zg  = 


C.  MODEL  CONSTRAINTS 

1.  Constraint  Formulation 

-  Ensure  weapon  availability  for  tasks  allocated  to  half-modules,  torpedo  tubes 
and  canister  capsules: 


ys,h,w\<h<m^ 


^b,wnZ  >0 


>0 


Z'V'V  <rj^^ 

/  j/  j  bwtpm  bw 


,heB'[,  P  m=\ 
weWf 


+ 


+ 
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-  Establish  the  allocation  at  the  bring  unit  level: 

yS  yS 

shwtpm  stpm  *> 


h-i  wGWi^\n^^>0 


\/s, t,p,m  s  ^  Sf  ,\<m<n^ 


ZyBT  Sr  yBC  _yB  \  t 

^hwtpm'^  2^  ^bwtpm~  hpm-’  \/ b  ,t ,  f  ,m\b  ^  ,l<m<p 


A.6.  + 

seSf  b<E^ 


\<  m<  n 


M 


Ensure  ah  (or  none)  of  the  task-parts  in  a  task  are  allocated: 


A7  'VE'^  -i-'VE'®  -t-VE 


B 

btpm 


seSf 


beB 


seSi  beB; 


^t,p,m\l<m<n^p  ,n,p,  >0, 

for  y  "G//"},m'  =  "l" 


Prevent  near-simultaneous  tasks  in  the  same  half-module: 


< 


A. 8.  ^  ^  ^  ^  !  ^shwtpm  ~ 

P  >n=l 

'  \w(PiVj 


ys,h,i 


\<h<m 


H 


Z  Z  ZZa 

P  m=l 


-BT  ^  BT 

w 

bwtpm  ’ 


\/b,i 


‘  WwJ 


Allocate  each  task-part  missile  requirement  or  take  a  penalty: 


A.10.  E  C  +  E  C+t^,..=i. 

seSj 


\ft,p,m 


\<m<n, 


M 


b^B 


Ensure  primary  and  back-up  task-part  relations: 


A.ii.  y:^+y;.^<x 


ses; 


\/s,t,m,m'^s  r 

for  p  =  "PR",p  ="BU" 
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A  17  F®  +F®  <1 


yb, t, m,m'\b  e  Bf  ,  \  <  m  <  ,  \  <  m'  <  nf, , 


for  p  =  "PR",p  ="BU" 


-  Ensure  primary  and  ready-spare  task-part  relations: 


M 


A  1 7  7®  =7® 

^stpm  ^stp'm^ 


A  1 4  =y^ 

r-v.A-r.  ^btpm  ^btp'm^ 


ys,t,m  s  e  ,\<m<n  , 


for  p  =  "PR'\p  ="RS" 


M 


beB,  ,\<m<pp. 


yb,t,m^ 

for  p  =  "PR",p'  =  "RS" 


-  Caleulate  firing  unit  engagement  in  targets  and  in  the  overall  alloeation: 


A.15. 

y  s,t,p,m,g\s  ,\<m<nfp 

A.16.  Z4<Zf, 

V5,g 

A.17.  7,;„<ZZ4, 

yb,t,p,m,g\beB^  ,\<m<nfp 

A.18.  ZZj<Zf, 

^b,g 

-  Limit  firing  units  for  restrieted  targets: 

A.19.  X^Zi+Z^^JS-r, 

s  b 

VgeG" 

-  Alloeate  all  (or  none)  of  the  tasks  in  restrieted  targets: 

A.20.  U,^^=U,^, 

'^t,t',g 

n,^>0,for  p^{'PR","C 
t,t'  eT^,t^t',g 

-  Use  of  expend  firing  units  for  primary  task-parts: 

< 

A.21.  X  Z  ZU.  for  p  =  "Pj;" 

seS'^^’’  tlsES!'  '«  =  l 


A.22.  X  Z 

bEBE”!’  t\bEBj  m=l 


for  p  =  ''PR'' 
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-  Ensure  the  allocation  of  all  missiles  for  non-separable  task-part  to  the  same  firing  unit: 


A.23.  £ 


m=l 


A-24.  £C=<"C. 


W=1 


\fs,t,p 


seSf,pelf\n:>2 


^b,t,pbeBf,pelf\n^>2 


-  Compute  absolute  deviation  from  average  number  of  residual  weapons  onboard  each 
firing  unit: 


A.25. 

h=l  w  i^seSJ  '«=! 

for  p  =  ”  PR " 

A.26. 

W  m=\ 

for  p  =  ”PR” 

A.27. 

^  ^  RN^  =  n^’^’^^RAVG 

A.28. 

RN^  - RAVG  > -RD^ , 

yseS\S^^P 

A.29. 

RN^  - RAVG  <  RD^ , 

^seS\S^^^ 

A.30. 

RN^  - RAVG  > -RDl^ , 

VbeB\B^^P 

A.31. 

RN^  - RAVG  <  RD^ , 

ybeB\B^^P 

-  Keep  track  of  areas  that  must  spread  primary  task-part  missiles: 

A.32. 

Z  ZC^raf, 

rseS/’  '"=1 

Vs  af  G  , 

for  p  =  ”PR” 

A.33. 

Z  ZC^/’^". 

'«=1 

Vb  G 

for  p  =  ”PR” 
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-  Keep  traek  of  areas  that  must  spread  baek-up  task-part  missiles: 


A.34. 


t\seS[  '«=1 


\/s 


a^eA^^,  for  p  =  "BU" 


A.35.  X 


tbeSj'  »«=1 


Vb 


for  p  =  "BU" 


Caleulate  residual  salvo  size: 


A.36.  ^  ^  ^ shwtpm-  ^shw-’ 


ieSf,  m=\ 


ys,h,w\<h<m^  ,  >  0,  for  p  =  "  PR 


A.37.  X  OL  ^  1. 


V^,/?  1  <  h  <  mf 


A.38.  n"-X  p  =  "PR' 


]beB,  ,  ^”=1 


\/b,w  >0,  for  p  =" PR 


A.39.  X 


BT 
b  ’ 


>0 


yb 


A.40.  nZ-X  Z  AiS,  S  O.t'';  P  =  "F3i",  Vi,w|«“>0,  for;,  =  "FF" 

JieS?’,  m=\ 


weW; 


-  TSC  speeifleations: 


A.41. 


V(s,  h,  w,  t,  p,  m)  eTSCf 


A.42.  XZp^=l, 


\/(b,  w,  t,  p,  m)  eTSCf^ 


A.43.  Xfi^=l, 


\/(b,  w,  t,  p,  m)  eTSCf^ 
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-  Submarine  back-up  allocation; 


A.44. 

A.45.  XZ,^=0, 

-  Force  levels  of  achievement  for  goals; 
A.46.  Z,<z:,(\  +  e,), 

AA7.  Z,>Z:.il-s,), 


b  i  ,w^W^, 

BT  /~v  1  _ '  _ '  M 

>  0,1  <  m  < 


for  p  =  "  BU", 


\/b,w,t,m 


b  i  B'^^  ,w^W^, 
>  0,1  <  m  <  n 


M 

tp 


for  p  =  "  BU", 


\/k'\k'<k, 
yk'\k'<k,  r,.=\ 


-  Binary  variables; 


A.48. 

\fs,h,w,t,p,ms  eS^  ,l<h<mf  ,weJ^  ,t  eT,r(^>0,peP,\<m<n^ 

A.49. 

yb,w,t,p,mbe^ ,weW^ ,n^  >0,t  eT,peP,\<m<n^ 

A.50. 

^b,w,t,p,mb&E(  ,w&Wl  ,n^>0,t&T,p&P,l<m<n^ 

A.51. 

\/s,g\seSf,geG 

A.52. 

zfs{o.i}. 

A.53. 

ZZ"  s{0,l}. 

yb,gbeBf,geG 

A.54. 

zfs{0,l}, 

ybeBj 
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A.55. 

V5,  t,  5  €  5f ,  t  €  r,  /?  €  7* 

A.56. 

'^b,t,p  b  eBf,teT,p  eP 

A.57. 

yseSf 

A.58. 

€{0,1}, 

ybeBj 

A.59. 

€{0,1}, 

yseSf 

A.60. 

€{0,1}, 

ybeBj 

-  Non-negative  Variables  (and  bounds,  if  applicable): 

A.61. 

0  <  7'^  <1 

^  —  ^stpm  — 

\/s,t,p,m s  eSj ,t  eT,p  eP,\<m<n^ 

A.62. 

yb,t,p,mbe^  ,t  eT,peP,l<m<n^ 

A.63. 

yt,p,m\^teT,peP,l<m<n^ 

A.64. 

\/s,h,w^s  €  Sf  ,\<  h  <  mf  ,we  , 

A.65. 

yb,w\b  e  B^  ,w  eW^  ,nZ  >0 

A.66. 

\/b,w\b^B^,w^W^,nZ>6 

A.61. 

TWf  >  0, 

yseS^ 

A.68. 

TWf  >  0, 

ybeBj 

A.69. 

IV 

o 

yseS^ 

A.70. 

IV 

o 

ybeBf 

A.71. 

RAVG>6 

A.72. 

>6 
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A.73. 


2.  Multi-Objective  Formulation 

We  consider  two  formulations  for  our  multi-objective  optimization  problem; 

1.  Hierarchical  formulation  (H);  For  eachA:e7f,  solve  the  following 
sequence  of  models; 


{H^) ;  Zl  =  min  (or  max) 

s.t.  (A.1)-(A.73) 


Remark;  To  be  more  precise,  we  write  all  these  models  as  the  following 
minimization  problems; 

iH,):  Z;=-/,  min-/,Z, 

s.t.  (A.1)-(A.73) 


2.  Combined-goal  (sometimes  referred  to  as  weighed-sum)  formulation  (C); 

(C);Z:=minZ„ 

s.t.  (A.1)-(A.45) 

(A.48)  -  (A.73) 


3.  Description  of  Formulation, 

•  Goals  #1  (Zl)  through  #8  (Zg)  seek  to  minimize  (or  maximize)  our  multiple 
objectives.  Zo  is  formulated  as  a  linear  combination  of  Zi,  ...,  Zg,  where  the  more 
important  the  goal,  the  higher  the  associated  weight. 

•  Constraints  (A.l)  to  (A.3)  ensure  that;  (a)  we  do  not  exceed  the  available  weapons  of 
each  type  in  any  half-module  on  a  ship,  and  (b)  we  do  not  exceed  the  available 
weapons  of  each  type  in  the  torpedo  room  and  in  the  canister  capsules  on  a 
submarine. 

•  Constraints  (A.4)  to  (A.6)  relate  the  half-module,  torpedo  room  and  canister  capsule 
variables  to  the  firing  unit  level.  In  particular,  (A.4)  tracks  the  ship  (if  any)  used  for  a 
specific  task-part  mission,  whereas  (A.5)  does  the  same  for  submarines.  (A.6)  ensures 
that  each  task-part  missile  requirement  is  allocated  once  (at  most). 
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•  Constraint  (A.7)  requires  the  alloeation  of  either  all  or  none  of  the  task-part  missiles 
in  a  task. 

•  Constraints  (A. 8)  and  (A.9)  avoid  alloeating  more  than  one  TLAM  per  half-module 
(for  ships)  and  more  than  the  number  of  torpedo  tubes  available  (for  submarines)  for 
task-parts  assoeiated  with  near-simultaneous  tasks. 

•  Constraint  (A.  10)  determines  whether  a  task-part  missile  requirement  has  been 
alloeated  or  not. 

•  Constraints  (A.  11)  and  (A.  12)  establish  the  relation  between  primary  and  back-up 
missions:  primary  and  back-up  task-parts  for  a  task  cannot  be  allocated  to  the  same 
firing  unit. 

•  Constraints  (A.  13)  and  (A.  14)  establish  the  relation  between  primary  and  ready-spare 
missions:  primary  and  ready-spare  task-parts  for  a  task  must  be  allocated  to  the  same 
firing  unit. 

•  Constraints  (A.  15)  and  (A.  17)  track  which  firing  units  have  been  used  in  the 
allocation  of  a  specific  target. 

•  Constraints  (A.  16)  and  (A.  18)  track  which  firing  units  have  been  used  in  the  overall 
allocation. 

•  Constraint  (A.  19)  limits  the  number  of  firing  units  assigned  to  a  restricted  target  not 
to  exceed  the  maximum  number  of  firing  units  allowed  for  that  target. 

•  Constraint  (A.20)  ensures  that  all  or  none  of  the  tasks  in  a  restricted  target  must  be 
allocated. 

•  Constraints  (A.21)  and  (A.22)  determine  if  an  expend  firing  unit  is  assigned  to  a 
primary  task-part  missile  requirement. 

•  Constraints  (A.23)  and  (A.24)  ensure  that  all  missile  requirements  for  a  non-separable 
task-part  are  allocated  to  the  same  firing  unit. 

•  Constraints  (A.25)  to  (A.31)  calculate  the  absolute  deviation  from  the  average 
number  of  residual  weapons  onboard  each  non-expend  firing  unit. 

•  Constraints  (A.32)  and  (A.33)  keep  track  of  the  number  of  firing  units  that  have  been 
allocated  primary  task-parts,  for  areas  so  designated. 

•  Constraints  (A.34)  and  (A.35)  keep  track  of  the  number  of  firing  units  that  have  been 
allocated  back-up  task-parts,  for  areas  so  designated. 

•  Constraints  (A.36)  to  (A.40)  calculate  the  residual  salvo-size.  In  particular,  (A.37) 
ensures  that  only  one  weapon  per  half-module  counts  for  computing  residual  salvo. 
(A.  3  9)  limits  the  maximum  count  for  weapons  in  the  torpedo  room  on  a  submarine  to 
the  number  of  torpedo  tubes. 

•  Constraints  (A.41)  to  (A.43)  satisfy  the  specifications  of  the  TSC. 
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•  Constraints  (A.44)  and  (A.45)  ensure  that  back-up  task-parts  are  not  allocated  to 
submarines  that  are  not  designated  for  back-up  tasking. 

•  Constraints  (A.46)  and  (A.47)  establish  levels  of  achievement  for  goals  that  are 
hierarchically  superior  to  the  goal  that  is  being  minimized  (or  maximized)  at  the 
incumbent  iteration.  (Only  for  the  hierarchical  formulation  (H).) 
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